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Abstract 

Motivated by recent evidence indicating that Quantum Einstein Gravity (QEG) 
might be nonperturbatively renormalizable, the exact renormalization group equa- 
tion of QEG is evaluated in a truncation of theory space which generalizes the 
Einstein-Hilbert truncation by the inclusion of a higher-derivative term (R 2 ). The 
beta-functions describing the renormalization group flow of the cosmological con- 
^ ■ stant, Newton's constant, and the i? 2 -coupling are computed explicitly. The fixed 

point properties of the 3-dimensional flow are investigated, and they are confronted 
with those of the 2-dimensional Einstein-Hilbert flow. The non-Gaussian fixed point 
predicted by the latter is found to generalize to a fixed point on the enlarged theory 
space. In order to test the reliability of the i? 2 -truncation near this fixed point we 
analyze the residual scheme dependence of various universal quantities; it turns out 
to be very weak. The two truncations are compared in detail, and their numerical 
predictions are found to agree with a suprisingly high precision. Due to the con- 
sistency of the results it appears increasingly unlikely that the non-Gaussian fixed 
point is an artifact of the truncation. If it is present in the exact theory QEG is 
probably nonperturbatively renormalizable and "asymptotically safe". We discuss 
how the conformal factor problem of Euclidean gravity manifests itself in the exact 
renormalization group approach and show that, in the i? 2 -truncation, the investi- 
gation of the fixed point is not afflicted with this problem. Also the Gaussian fixed 
point of the Einstein-Hilbert truncation is analyzed; it turns out that it does not 
generalize to a corresponding fixed point on the enlarged theory space. 
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1 Introduction 



Recently a lot of work on Quantum Einstein Gravity (QEG) went into constructing an 
appropriate exact renormalization group (RG) equation |], ||, finding approximate solu- 
tions to it [0, |], ||, ||, [7j , and exploring their implications for black hole physics |8], |9| and 
cosmology [|l(J. In particular, strong indications were found that QEG might be nonper- 
turbatively renormalizable. If so, it could have the status of a fundamental, microscopic 
quantum theory of gravity. 

The basic tool used in these investigations is the effective average action and its exact 



RG equation (TT| . It is a continuum analog of Wilson's lattice renormalization group of 
iterated block spin transformations Both in quantum field theory and statistical 

mechanics the idea is to integrate out all fluctuation modes which have momenta larger 
than a certain infrared (IR) cutoff k ("fast degrees of freedom"), and to take account 
of those modes in an implicit way by the modified dynamics which they induce for the 
remaining fluctuations with momenta smaller than k (the "slow degrees of freedom"). 
In field theory this "renormalized" dynamics is encoded in a scale dependent effective 
action, T^, whose dependence on the cutoff scale k is governed by a functional differential 
equation referred to as the "exact RG equation" [O. This equation gives rise to a flow 
on the space of all actions ("theory space"). The functional Tf, defines an effective field 
theory valid near the scale k; evaluated at tree level, it describes all loop effects due to 
the high-momentum modes. The effective average action can be thought of as a kind of 
microscope with a variable resolution. At large k, the physics at short distances £ = 1/k 
can be read off directly from 1^; at small k we see a coarse-grained picture suitable for a 
simple description of structures with a large characteristic length scale i = 1/ k . 

The effective average action T^, regarded as a function of k, interpolates between the 
ordinary effective action T = lim^o and the bare (classical) action S which is ap- 
proached for k — > oo. The construction of begins by adding a IR cutoff term A^S 
to the classical action entering the standard Euclidean functional integral for the gener- 
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ating functional W of the connected Green's functions. The new piece A k S introduces 
a momentum dependent (mass) 2 -term lZ k {p 2 ) for each mode of the quantum field with 
momentum p. For p 2 ^> k 2 , the cutoff function TZ/,(p 2 ) is assumed to vanish so that 
the high- momentum modes get integrated out unsuppressed. For p 2 <C k 2 , it behaves 
as TZk{p 2 ) oc k 2 ; hence the small- momentum modes are suppressed in the path integral 
by a mass term oc k 2 . Apart from a correction term which is known explicitly ||11|| , the 
effective average action is given by the Legendre transform of the modified generating 
functional W k . 

From this definition one can derive the exact RG equation obeyed by T^. In a slightly 
symbolic notation it is of the form 

1 



k d k Y k = - Tr 



(ri 2) + ^ fe (-A)) 1 kd k K k (-A) 



The RHS of this equation is a kind of "beta-functional" which summarizes the beta- 
functions for infinitely many running couplings. Geometrically, it defines a vector field on 
theory space, the corresponding flow lines being the RG trajectories k i— > r&. 

The functional T k enters this vector field via its Hessian T k 2 \ i.e. the infinite-dimensional 
matrix of all second functional derivatives of T k with respect to the dynamical, i.e., non- 
background fields. 

In eq. Ql.lQ the c-number argument of TZ k is replaced with the operator —A. The 
discrimination of high- "momentum" vs. low- "momentum" modes is performed according 
to the spectrum of this operator, i.e. p 2 is an eigenvalue of —A. In simple theories where 
no gauge or diffeomorphism invariance needs to be respected, A is the free Laplacian, 
A = d^d^, whose eigenmodes are momentum eigenstates in the usual sense of the word. In 
Yang-Mills theory |TJ|, [TJ| it has proven convenient to use the background field formalism 



lq , 17 and to set A = DnD^ where is the covariant derivative in the background 



field. The background field technique plays a dual role in this context. Using a background 
gauge fixing term makes T k a gauge invariant functional of its argument, and using D M in 
the cutoff leads to a flow equation of the relatively simple type ( |1.1|) , similar to non-gauge 
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theories. (The RG equation resulting from A = D^D^ with constructed from the 
dynamical gauge field is quite unwieldy.) 

Along similar lines an effective average action for o?-dimensional Euclidean QEG has 
been constructed in ref . Jl| , and the corresponding flow equation has been derived. Leav- 
ing the Faddeev-Popov ghosts aside for a moment, the gravitational effective average ac- 
tion, r^g^, g^v], is a functional of two different metrics, the "ordinary" dynamical metric 
g^u and the background metric g^. The usual effective action T^j,] is recovered by tak- 
ing the limit k — > of the functional T^g^} = Tk[g^i,, g^ = g^u] in which the two metrics 
are taken equal. Thanks to the background gauge fixing condition this construction leads 
to a functional r^g^] which is invariant under general coordinate transformations. 

One of the many advantages which the exact RG approach has in comparison to 
the standard canonical or path integral quantization is that it offers a very natural and 
intuitive nonperturbative approximation scheme. By truncating the theory space one can 
obtain approximate solutions to the RG equation which do not need a small expansion 
parameter. The idea is to project the RG flow from the "huge" infinite-dimensional space 
of all actions onto some smaller, typically finite dimensional subspace which is easier to 
handle. In this way the functional RG equation for becomes a system of ordinary 
differential equations for a (finite) set of coupling constants which have the geometrical 
interpretation of coordinates on the subspace. It is clear that in applying this strategy 
the key problem is finding the "relevant" subspace which contains the essential physics. 

In a first attempt at solving the gravitational RG equation jy] the flow has been pro- 
jected onto the 2-dimensional subspace of theory space which is spanned by the invariants 
J d d x sjg and J d d x ^/g~R. This is the so-called Einstein-Hilbert truncation defined by the 
ansatz 

Fk[g, g) = (W^Gk) 1 J d d x ^fg {—R(g) + 2Afc} + classical gauge fixing. (1.2) 

The two running couplings involved are the running Newton constant and the running 
cosmological constant The similarity of ( |1 . 2|) to the action of classical General Rela- 
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tivity is accidental in a sense; improved truncations would include both higher powers of 
the curvature and nonlocal terms [[TJ|, [J. In ( |1.2| ) only the gauge fixing term depends on 
g^ v ; it vanishes when we set g^ v = g^ v . 

The scale dependence of Gk and A& is most conveniently visualized as a flow in the 
\-g— plane where gk = k d ~ 2 Gk and A& = Xk/k 2 are the dimensionless Newton constant 
and cosmo logical constant, respectively. Using the original cutoff of "type A" |j] the 
system of equations for gk and A& was derived in |J and solved numerically in ||. In 
0] a new cutoff of "type B" was introduced and the corresponding flow equations in the 
Einstein-Hilbert truncation were derived. (The "type B" cutoff is convenient if one uses 
the TT-decomposition of the metric The fixed point properties of these equations 

were first discussed in |20| and ||, and analyzed in detail in || and ||. One finds that the 
RG flow in the X-g— plane is governed by two fixed points (A*, g*)\ a trivial or "Gaussian" 
fixed point at (A*, g*) = (0, 0), and a non-Gaussian fixed point with A* ^ and g* ^ 0. 

In order to appreciate the importance of the non-Gaussian fixed point we recall what 
it means to "quantize" a theory in the average action approach. One picks a bare action 
S and imposes the initial condition IV = S at the ultraviolet (UV) cutoff scale k, uses the 
RG equation to find Tk at all lower scales k < k, and finally sends k — > and k — > oo. 
A fundamental theory has the property that the "continuum" limit k —>■ oo actually 
exists after redefining only finitely many parameters in the action. This is the case in 
perturbatively renormalizable theories , but there are also examples of perturbatively 
nonrenormalizable theories which possess a limit k — > oo The continuum limit of 

those "nonperturbatively renormalizable" theories is taken at a non-Gaussian fixed point, 
i.e. the theory is defined by the set of RG trajectories which leave the fixed point when we 
lower k. These trajectories span the UV critical hypersurface of the fixed point, <Suv- If 
it is finite dimensional, the quantum theory thus constructed has only finitely many free 
parameters and therefore keeps its predictive power even at arbitrarily large momentum 
scales. This behavior is to be contrasted with an effective field theory which, at high 
energies, typically contains an increasing number of free parameters which must be taken 
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from the experiment. 

In his "asymptotic safety" scenario Weinberg |2B|, 2J] conjectured that a fundamen- 



tal quantum field theory of gravity could perhaps be constructed nonperturbatively by 



taking the continuum limit at a non-Gaussian fixed point P5 |. While originally this idea 
could be implemented in d = 2 + e dimensions only, the recent results coming from the 
effective average action strongly support the hypothesis that this fixed point exists also 
in 4 dimensions. Within the Einstein-Hilbert truncation, the existence of a suitable non- 
Gaussian fixed point is definitely established by now; the crucial question is whether it 
is the projection of a fixed point present in the exact theory or merely an artifact of the 
approximation. 

Let us assume for a moment that the fixed point indeed exists in the exact 4-dimen- 
sional theory and that we define QEG by taking the k — ► oo limit there. Then, since 
(Afc, gk, ■ ■ ■) approaches (A*, g*, ■ - •) for k — > oo, the dimensionful couplings behave as 

G k ^gjk\ \ k ^Kk 2 , •■■ (1.3) 

for large k. Obviously G%. vanishes for k — > oo. At least as far as this coupling is con- 
cerned QEG is asymptotically free similar to Yang-Mills theory. The l/fc 2 -dependence 
of the running Newton constant will lead to a characteristic momentum dependence of 
the cross sections for graviton-graviton scattering and graviton mediated matter-matter 
scattering. Because of this characteristic momentum dependence, QEG could be distin- 
guished experimentally from alternative theories of quantum gravity such as string theory, 
at least in principle. 

Let us generalize the standard definition of the Planck mass, mp\ = G^ 1 ^ 2 , and intro- 
duce the running Planck mass 

M m (k) = l/y/(h. (1.4) 

At the laboratory scale Mp\(k) reduces to mpi, most probably, and its dependence on k 
is negligible. However, in the fixed point regime k — >• oo, the asymptotic freedom of G k 
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implies that Mp\(k) is proportional to the scale k itself: Mp\(k) = kj-Jg~l. This shows 
that the running Planck mass is a rather elusive "barrier" which never can be jumped 
across in any experiment. If we analyze a system with a probe of increasing momentum 
k we will always push the running Planck mass ahead of us and never reach it. 

Also the standard constant Planck mass, defined more precisely in terms of the IR 
value of Gk, 

mpi = [G fe=0 ]- 1/2 , (1.5) 

plays an important role in QEG, similar to that of Aqcd in QCD. According to the 
numerical solutions of the A- g-system || , rap\ marks the lower boundary of the asymptotic 
scaling region. Near k = mp\ there is a crossover from the scaling laws (|1.3|) of the non- 
Gaussian fixed point to those of the Gaussian fixed point 

According to the UV scaling laws (|1.3|) the dimensionful cosmological constant diverges 
for k — > 00 proportional to k 2 . This has an interesting geometrical interpretation. Let us 
consider the ^-dependent, effective field equations implied by the truncation ansatz ( |1.2| ) 
with g^y = g^ v for d = 4. They happen to coincide with the familiar vacuum Einstein 
equations with the cosmological constant replaced by the scale-dependent quantity X^. 

Rfus - ^9fiu R = -Afe g^u . (1.6) 

Since A& is the only quantity which sets a scale, every solution to (|1.6|) has a typical radius 

of curvature r c (k) oc 1/ vAfc- (F° r instance, the maximally symmetric S" 4 -solution has the 

radius r c = r = \/3/Afc.) The ^-dependence of the solutions and in particular of r c should 

be interpreted as follows. If we want to explore the spacetime structure at a fixed length 

scale I = 1/k it is most convenient to use the action r^g^] at k = \jl because for 

this, and only this, functional a tree level analysis is sufficient to describe the essential 

physics at this scale, including all quantum effects. Hence, when we observe spacetime 

with a "microscope" of resolution £, we will see an average radius of curvature given by 
1 A similar crossover was already known to occur in Liouville quantum gravity pq|. 
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r c (£) = r c {k = l/£). Once £ is smaller than the (standard) Planck length l?\ = m p ^ we 
are in the fixed point regime ( |1.3| ) so that r c {k) oc 1/k, or 

r c (£)^£ (1.7) 

Thus, when we look at the structure of spacetime with a microscope of resolution £, the 
average radius of curvature which we measure is proportional to the resolution itself. If we 
want to probe finer details and decrease £ we automatically decrease r c and hence mcrease 
the average curvature. Spacetime seems to be more strongly curved at small distances 
than at larger ones. The scale-free relation ( |1.7| ) suggests that at distances below the 
Planck length quantum spacetime is a kind of fractal with a self-similar structure. It has 
no intrinsic scale. 

Before we continue a remark might be in order on which theory precisely we refer 
to as "Quantum Einstein Gravity" or "QEG". While flow equations can also be used 
in the effective field theory approach to quantum gravity |27], |IJ, in the present context 
"QEG" stands for the fundamental theory whose continuum limit k — > oo is taken at the 
non-Gaussian fixed point. This theory has dim(iSuv) free parameters, and fixing these 
parameters amounts to picking a specific trajectory k i— > in the full theory space. For 
k — > oo this trajectory hits the fixed point action T*, regarded as the collection of its 
infinitely many dimensionless coordinates on theory space. The fixed point action r* 
corresponds to the "bare" or "classical" action in conventional field theory. However, 
unlike the latter r* is not put in by hand but is rather derived with the help of the 
RG equation. The usual canonical or path integral quantization is always based upon 
a "prejudice" about what the classical action is. In the asymptotic safety scenario the 
"classical" action is fixed instead by the condition of nonperturbative renormalizability; 
it cannot be guessed by simple power-counting, symmetry, or invariance arguments, but 
the effective average action provides a computational framework to determine it. 

The really crucial property which defines QEG is the requirement of diffeomorphism 
invariance. Before we can write down a flow equation we must declare what the theory 
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space is on which the renormalization group is supposed to operate. In the case of QEG it 
is defined to be the space of functional r[^] depending on a nondegenerate, symmetric 
rank-2 tensor field in a diffeomorphism invariant way. Leaving technical details aside, the 
theory space then fixes the flow equation which in turn determines the RG trajectories 
and the fixed points. 

So far our discussion referred to the exact theory on the full theory space. If we 
project on the subspace spanned by the Einstein-Hilbert truncation it is clear that also 
the fixed point action must be of the Einstein-Hilbert type. However, we emphasize 
that this is a trivial consequence of the simple truncation we have chosen, and we have no 
reason to believe that the exact is of the Einstein-Hilbert type, too. In fact, within the 
more complicated truncation of the present paper r* receives corrections which go beyond 
the Einstein-Hilbert form. Hence "Quantum Einstein Gravity" does not mean that the 
Einstein-Hilbert action is the bare action to be quantized. In this respect our approach is 
different from canonical quantum gravity, along the lines of Ashtekar's program |E3|, for 
instance. (It is intriguing that also in this context remarkable finiteness properties have 
been proven recently fl29fl .) 

Clearly it is a highly attractive idea that there could be a nonperturbatively renormal- 
izable field theory of the metric field so that there is no longer any conceptual need for 
leaving the framework of quantum field theory in order to arrive at a consistent micro- 
scopic theory of quantum gravity. Therefore every effort should be made to show that the 
non-Gaussian fixed point found in the Einstein-Hilbert truncation is not just an artifact 
of this approximation. 

In ref. we therefore started an extensive analysis of the reliability of the Einstein- 
Hilbert truncation near the fixed point. There the strategy was to use the scheme depen- 
dence of universal quantities in order to get a first idea about the precision which can be 
achieved with this truncation. Here "scheme dependence" refers to the dependence on the 
details of the cutoff procedure, i.e. on the shape of the function H k {p 2 ). By definition, 
universal quantities are exactly scheme independent in the exact theory, but they might 
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acquire some scheme dependence once we make approximations. The level of this residual 
scheme dependence can serve as a measure for the quality of the approximation. Typi- 
cal universal quantities are the critical exponents of fixed points and, as we argued, the 
product g* A*. The upshot of our analysis was that the Einstein-Hilbert truncation seems 
to provide a description that is much more reliable and precise than originally hoped for, 
and that it would be very hard to understand the approximate scheme independence we 
found if the fixed point was just due to a misleading approximation. 

These results are certainly very encouraging, but it is clear that the ultimate justifi- 
cation of a truncation ansatz consists of adding further terms to it and verifying that its 
predictions do not change much. In the present paper we take a first step in this direction 
and add one further invariant constructed from to the ansatz. 

Which invariant should we take? In standard renormalized perturbation theory where 
(at least implicitly) the k — > oo limit is taken at the Gaussian fixed point, the relative 
importance or "relevance" of the various field monomials is measured by their scaling 
dimensions at the Gaussian fixed point, i.e. by their canonical dimensions simply. Since 
at the non-Gaussian fixed point the anomalous dimensions are large we have no similarly 
simple guide line at our disposal, and a priori all invariants are equally plausible. To get 
a first idea about what happens away from the Einstein-Hilbert subspace we shall include 
the higher-derivative invariant J d d x ^fg~R 2 and study the RG flow in the ".R 2 -truncation"0 

Fk[g,g] = J d d x yfg {(167rG fc ) _1 [-R{g)+2X k ] + p k R 2 (g)} + classical gauge fixing(1.8) 

Its truncation subspace is 3-dimensional, with coordinates G, A and the new coupling j3. 

It is well known that beyond J d d x ^fg~R 2 there exist twc0 more (curvature) 2 -invariants: 

J d d x ^fgR iiV R >iV and J d d x ^fgR VLVpa R iiVpa '. In a standard perturbative calculation near 

the Gaussian fixed point consistency would require to include them along with the i? 2 -term 

because they all have the same canonical dimension. As for the non-Gaussian fixed point, 

we have no a priori information from general principles about the relative importance 
2 Our conventions are R° ' piiv = -d v T^ p H , R pv = R a pai/ , R = g^R^v 

3 Except in d — 4 where one invariant can be eliminated by virtue of the Gauss-Bonnet identity. 
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of the three terms. Since anyhow the best we can do is to take a "step into the dark", 
without knowing whether we walk in the most "relevant" direction, we shall omit the 
other two invariants here. Including them would go far beyond the present calculational 
possibilities, in particular since it would require a much more complicated projection 
technique 0. 

In this paper we shall derive the (extremely complicated) 3-dimensional RG equations 
of the -R 2 -truncation, and we shall use them in order to investigate the fixed points of the 
flow. Our main results will be the following: a) The Gaussian fixed point of the Einstein- 
Hilbert truncation does not generalize to a fixed point of the i? 2 -truncation. b) The non- 
Gaussian fixed point does indeed generalize to a fixed point of the _R 2 -truncation, and the 
X-g— projection of this fixed point is described almost perfectly by the Einstein-Hilbert 
truncation. Within the (weak) residual scheme dependence, the fixed point properties 
are almost insensitive to the inclusion of the i? 2 -invariant. This is further strong evi- 
dence against the theoretical possibility that the non-Gaussian fixed point is a truncation 
artifact. 

In the second part of the paper we shall address a very important general problem 
which is of a more technical nature. It is related to a notorious disease of standard 
Euclidean quantum gravity: the conformal factor problem. In setting up the truncated 
RG equation the cutoff function IZk (actually a matrix in field space) is adapted to the 
truncation in such a way that, for p 2 <C k 2 , the inverse propagator of every massless 
mode, p 2 , is replaced by p 2 + k 2 . A problem arises if there are modes, such as those 
of the conformal factor in the Einstein-Hilbert truncation, which have a negative kinetic 
energy, i.e. their inverse propagator is — p 2 . In |T[ it has been argued that for these modes 
also the sign of TZ/. should be reversed so as to obtain the regularized inverse propagator 
— (p 2 + k 2 ). While there is little doubt that this procedure is correct for the Einstein- 
Hilbert truncation, it leads to the seemingly paradoxical situation that in the Euclidean 
path integral the modes of the conformal factor are enhanced rather than suppressed in 
the IR 0. 
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Contrary to the Einstein-Hilbert truncation, the i? 2 -truncation yields a functional 
rfe^i/] which is bounded below and, as we shall see, gives positive kinetic energy to all 
modes, provided one stays close to the UV fixed point. Hence our investigation of the non- 
Gaussian fixed point is not plagued by the conformal factor problem, and the construction 
of the cutoff becomes straightforward. This advantage is an independent motivation for 
studying the i? 2 -truncation. 

In the usual perturbative approach higher-derivative theories of (Lorentzian) gravity 
are notoriously problematic as far as causality and unitarity are concerned J3U], ^T|. While 
in d = 4 the most general (curvature) 2 -theory, when expanded about flat space, is power- 
counting renormalizable it suffers from excitations with, classically, negative linearized 
energy and, quantum mechanically, a wrong-sign residue of the propagator leading to 



a state space containing negative- norm states [p0| . These "ghosts" have masses of the 
order of the Planck mass. Correspondingly, if a truncation ansatz is of the (curvature) 2 - 
type the fc-dependent effective propagator ^Tk/Sgdg)^ 1 = (T k evaluated for flat 
space, has similar ghosts with masses oc Mp\(k). By itself this does not indicate any real 
problem because generically flat space anyhow is not a solution of the effective equation 
of motion. Compared to the perturbative quantization the potential problem of ghost 
excitations manifests itself in the Euclidean RG approach in a conceptually different, more 
tractable manner. Here the linearization is performed about the backgrounds needed for 
the projection procedure, not about flat space. For a well-defined computation of the RG 
trajectories on a certain fc-interval it is sufficient that the (truncated) Tk gives positive 
linearized action to all modes contributing to the RG running in this interval. In our 
calculation this will indeed be the case for k large enough where the truncation is believed 
to be reliable ||. The much more subtle issues related to a Lorentzian interpretation of 
the theory and its causality properties can be addressed only once a complete trajectory, 
valid down to the IR, and in particular the precise form of the fixed point action is 
known. From all what we can tell now the exact QEG could very well be "causal" in an 
appropriate sense. 
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A brief summary of some of the results derived in the present paper appeared in ||, 
and an informal introduction to the older work can be found in |52| . In the present paper 
we focus on pure gravity. The gravitational average action with matter fields included 
was discussed in JF| and |33| . The gauge fixing dependence of the original formulation [1] 
was investigated in [33] and [35]. An incomplete higher-derivative calculation was begun 



in where the running of the i? 2 -couplings was neglected, however, and no conclusions 
about the fixed point could be drawn. 

The remaining sections of this paper are organized as follows. In Section [2] we review 
some general properties of the exact RG equation which will be needed later on. Section 
[5] is devoted to the construction of cutoffs which are adapted to a specific truncation; in 
particular the complications due to the conformal factor problem will be discussed there. 
In Section |] we derive the system of RG equations which results from the _R 2 -truncation, 
and in Section |5| we analyze its fixed point structure. We discuss the fate of the Gaussian 
fixed point which is present only in the Einstein-Hilbert truncation, and we reanalyze the 
non-Gaussian fixed point in the more general setting. In Section |] the positivity properties 
of the truncated action functional, its Hessian, and the cutoff operator are investigated; 
in particular we show that our analysis of the non-Gaussian fixed point is not affected 
by the conformal factor problem. The conclusions are contained in Section [7[ Many 
important technical results, including the coefficients occuring in the rather complicated 
beta-functions of the \-g-f3— system, are tabulated in various appendices. 



2 The exact RG equation 

In this section we briefly review the construction of the "type B" RG equation for quantum 
gravity performed in ref. [0] to which we refer for the details. We start from a scale- 
dependent modification of the generating functional for the connected Green's functions, 
Wk- It is defined by the following Euclidean functional integral: 

exp {W4 [sources]} = / Vh^ VC* 1 VC^ exp [S[g + h] 
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-S gi [h; g] - S gh [h, C, C- g] - A k S[h, C, C; g] - S source ] . (2.1) 

In ( |2.1| ) we use the background gauge fixing technique which necessitates the decomposi- 
tion of the full quantum metric 7^ into a fixed background metric g^ and a fluctuation 
variable h^ u : 7^ (x) = g^(x) + h^ u (x). It allows us to replace the integration over 7^ by 
an integration over h^ v . Furthermore, and C 1 are the Faddeev-Popov ghosts of the 
gravitational field. 

The first term of the action, S[g + h], is the classical part, which is assumed to be 
invariant under general coordinate transformations. For the time being, we also assume 
that it is positive definite, S > 0. The gauge fixing term is given by 

Sgf[M = ^f d d x^gg^F,[gMFA9A • (2.2) 

In the present paper we use the linear gauge condition F^g.h] = \/2k T^\g] h a p with 
T^\g] = 5^g ai D 1 — -g^D^, which amounts to a background version of the harmonic 
coordinate condition. Here we introduced the constant k = (32nG)~^ where G is the 
bare Newton constant. Moreover, denotes the covariant derivative constructed from 
the background metric g^, while we shall write for the covariant derivative involving 
the complete metric 7^. is the Faddeev-Popov ghost action resulting from the above 
gauge fixing. 

Furthermore, A^S and 5 sour ce are the cutoff and the source action, respectively. A k S 
provides an appropriate infrared cutoff for the integration variables and will be discussed 
in detail in a moment; S^m-ce introduces sources for the fields h^, C M and C^. 

Next we decompose the gravitational field according to (see e.g. |19|] ) 

V = h lu + D£ u + + D^Dyd - -g^D 2 a + X -g^ • (2-3) 

In order to obtain this "TT-decomposition" one starts by writing h^ u as a sum of its 
orthogonal parts: h^ v = + h^ u + h^. Here h/E v , hh v and hj£ represent the transverse 
traceless, longitudinal traceless and pure trace part, respectively. Introducing two scalar 
fields and a, and a transverse vector field the tensors UFj v and h^ u can be expressed by 
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tfL = Q^A/d and h%, = hff+hfi with hf v = and h% = D^a-g^a/d. 

Thereby we end up with eq. fl2.3|) . In the following the components of h^ u thus introduced 
will be referred to as the "component fields" . They obey the relations 



$"7& = 0, D^ u = 0, D% = 0, 4> = 9,uh^. (2.4) 

Obviously the complete field h^ v receives no contribution from those £ M - and a-modes 
which satisfy the Killing equation 

D„t + A,& = (2.5) 

and the scalar equation 

D li D v B-±g lu ,D 2 a = 0, (2.6) 

respectively. Such "unphysical" £ M - and a-modes have to be excluded from the functional 
integral and all subsequent calculations 0. Having a closer look at the scalar equation 
( O ), one recognizes that there is a one-to-one correspondence between the nonconstant 
solutions of (|2.6|) and the purely longitudinal, or proper, conformal Killing vectors (PCKV) 
C M . They are related via C M = D^d. 

Likewise we decompose the ghost and the antighost into their orthogonal components: 

C li = C% + D ll fj, C» = C T » + D»rj. (2.7) 

Here Cj and C Tm are the transverse components of and C M : -D M C*J = 0, D m C Tm = 
0. Furthermore, the scalars f) and ff parametrize the longitudinal part of and C M , 
respectively. The constant fj- and rj-modes represent unphysical modes which have to be 
excluded. 

For calculational convenience we now introduce new variables {£ fl ,&,f],T)} replacing 
a, f), rf}, by means of the momentum dependent (nonlocal) redefinitions 

f " = V-D 2 -Rlc" C 



-D 2 r) , r] = \/ -D 2 jf. (2.8) 
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Here the operator Ric maps vectors onto vectors according to (Ric f) M = R^Vy. In accor- 
dance with the decompositions ( 2T3 ), ( |2.7| ) and the redefinitions Q2.8| ) we then perform the 



combined transformation of integration variables h^ u — > {h? v , £ M , a, 0}, C7 M — >■ {Cj, fj}, 
C^ L — ► {C* Tm , 77} in the functional integral (|2.1|) . The Jacobian induced by this change of 
variables is such that it boils down to an umimportant constant if Einstein backgrounds, 
characterized by R^ u = Cg^ u with C a constant, are inserted into (|2.1| ). (For the general 
case see ref. @.) 

Let us now come to the "type B" cutoff term A k S. At the component field level, it is 
a sum of inner products, 

A k S[h,C,C;g]=\ £ (K k ) Ct(2 ( 2 ) + ~ £ (72*)** V*) (2-9) 



Ci>C2e/i Vi iV'a s/2 

with the index sets Ji = {/i T , £, a, 0}, J 2 = {C T , C T , Vi v}- At this stage of the discussion 
it is not necessary to specify the explicit structure of the cutoff operators lZ k acting on the 
component fields. In order to provide the desired suppression of low-momentum modes, 
these operators must vanish for p 2 /k 2 — > 00 (in particular for k — > 0) and must behave 
as TZk —> Z k k 2 for p 2 /k 2 —>■ 0. (The meaning of the constant Z k will be explained later.) 
Furthermore, they have to satisfy certain hermiticity conditions ||. 

Now we are in a position to construct the effective average action r*.. It is defined as 
the difference between the Legendre transform of W k at fixed g^ v , denoted T k \h,v,v;g], 



and the cutoff action with the classical fields inserted 37, 14 



T k [g, g, v, v] = T k [g- g, v, v; g] - A k S[g - g, v, v; g] . (2.10) 

Here the classical fields represent the (fc-dependent) expectation values of the quantum 
fluctuations: h^ u = (h^), = (C^), t> M = (C 1 '). They are obtained in the usual way 
as functional derivatives of W k with respect to the sources. In ( |2.10| ) we expressed h^ u 
in terms of the classical counterpart g^ v of the quantum metric 7^ = g^ v + h^ u which, 
by definition, is given by g^ v = g^ iu + h^. The classical analogs of the components 
{h T , £, a, 4>, C T , C T , fj, rf) will be denoted (h T , £, a, </>, v T , v T , g, g). 
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The exact RG equation for the effective average action describes the change of 
induced by an infinitesimal change of the scale k. Introducing the RG "time" t = In k, it 
can be derived from the t-derivative of the functional integral fl2.1|) . It takes the form 



d t T k [g,g,v,v] 



-Tr' 



E ( r ? } & 9, «, v] + K k ) ^ dt (K k \ 2(i 



Cl,C2€/l 



+-Tr' 
2 



£ ( r i 2) t^, 



v,v] +TZ 



-1 



/,. j d t i'Ri,.) 



^1,-026^2 



. f2.ll) 



(2) 

Here T k denotes the Hessian of T k with respect to the component fields. Furthermore, 
we wrote (TZk)^ = (^-fc)(Ci)(C2>) (T^k)^^ = {T^k)^)^) and introduced the index sets 
Ii = {h T ,^,a,(j)} , li = {v T , v T , g, g}. Furthermore, the primes at the traces indicate 
that all unphysical and cr-modes, characterized by eqs. ( |2~5| ) and (p.6|), are to be 
excluded from the calculation of the traces. 



3 Truncations and their adapted cutoffs 
3.1 Truncating the ghost sector 

In concrete applications of the exact RG equation one encounters the problem of dealing 
with an infinite system of coupled differential equations. Usually it is impossible to 
find an exact solution so that we are forced to rely upon approximations. A powerful 
nonperturbative approximation scheme is the truncation of theory space, which means 
that only a finite number of couplings is considered and the RG flow is projected onto 
a finite-dimensional subspace of theory space. In practice one proceeds as follows. One 
makes an ansatz for that comprises only a few couplings and inserts it on both sides of 
eq. ( |2.11|) . By projecting the RHS of this equation onto the space of operators appearing 



on the LHS one obtains a finite set of coupled differential equations for the couplings 
taken into account. 
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Given an arbitrary truncation it is not clear a priori whether it is sensible and leads to 
at least approximately correct results. In this respect the modified BRS Ward identities 
satisfied by the exact T k [0 are of special importance, since only those truncations which 
are (approximately) consistent with them can be reliable. In [|l| it was shown that under 
certain conditions truncations of the form 



T k \g,g,v,v) = T k [g] + T k [g,g] + S g{ [g -g;g] + S gh [g - g,v,v; g] 



(3.1) 



which neglect the RG running in the ghost sector are approximate solutions to the Ward 
identities for the exact T k . Here T k [g] is defined as 



r k \g] = T k \g,g, 0,0] 



(3.2) 



and T k [g,g] encodes the quantum corrections of the gauge fixing term. (For the details 
we refer to fl], §.) Inserting the ansatz ( j3~Tl) into the exact evolution equation (|2.11|) 
leads to a truncated RG equation which describes the RG flow of T k in the subspace of 
action functionals spanned by (|3.1| ). The equation governing the evolution of the purely 
gravitational action 



Fk[g,g) = T k [g,g, 0,0] = T k [g] + S gi [g - g; g] + T k [g,g] 



(3.3) 



takes the form 



d t T k [g,g] = -Tr' 



_ T v * C1C2 



Cl,C2 6/i 



1 . 

+ 2 Tl 



Here and are the Hessians of T k [g,g] and S g h[h, v, v;g] with respect to the grav- 
itational and the ghost component fields, respectively. They are taken at fixed g^. 



y(2) 



(3.4) 
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3.2 Construction of the cutoff, and the conformal factor prob- 
lem 

In order to obtain a tractable evolution equation for a given truncation it is necessary 
to use a cutoff which is adapted to this truncation but still has the desired suppression 
properties for a class of backgrounds which is as large as possible. 

A convenient cutoff which is adapted to the truncation ansatz can be found in the 
following way [[l], |7j. Given a truncation, we assume that for g = g the kinetic operators 
of all modes with a definite helicity can be brought to the form (T k )ij = fij(—D 2 , k, . . .) 
where {fij} is a set of c-number functions and the indices i,j refer to the different types 
of fields. Then we choose the cutoff operator 7Z k in such a way that the structure 

(if > + n k )ij = fa {-D 2 + k 2 R(°\-D 2 /k 2 ), k,...) (3.5) 

is achieved. Here the function R^(y), y = —D 2 /k 2 , describes the details of the mode sup- 
pression; it is required to satisfy the boundary conditions RS '(0) = 1 and lim R^(y) = 0, 
but is arbitrary otherwise. By virtue of eq. ( |3.5|) , the inverse propagator of a massless field 
mode with covariant momentum square p 2 = —D 2 is proportional to p 2 + k 2 R^(p 2 /k 2 ) 
which equals p 2 for p 2 ^> k 2 and p 2 + k 2 for p 2 <^k 2 . This means that the small-p 2 modes, 
and only those, have acquired a mass oc k which leads to the desired suppression. 

In order to see the potential problems of the rule ( |3.5| ) let us be more specific and 
assume that the functions fij are linear in D 2 and contain no constant term. Then, after 
diagonalizing with respect to the field indices, 1^ decomposes into a set of (massless, 
by assumption) inverse propagators z k p 2 with a running wave function normalization z k . 
In this diagonal basis TZ k is diagonal, too. It is of the form 

n k {p 2 ) = Z k k 2 R^{^J . (3.6) 

A priori Z k is a free constant, but when we apply the rule (|3.5| ) we are forced to set 
Z k = z k for each mode. Only then the propagator and the cutoff combine in the right 
way, leading to the modified inverse propagator z k [p 2 + k 2 R^\p 2 / k 2 )\. 



18 



The choice Z k = Zk is certainly the correct one if z& is positive. This is indeed 
the case in the familiar unitary theories on flat spacetime. In QEG, however, there are 
truncations, the Einstein-Hilbert truncation, for instance, which give a negative kinetic 
energy to certain modes (p of the metric. In particular, in the Einstein-Hilbert truncation, 
the conformal factor (p = <ft has z k < 0. 

The important question is how Z%. should be chosen when z k is negative. If we continue 
to use Zk = Zk, the RG equation is still well-defined because the inverse propagator 
— |zfc|[p 2 + k 2 R(°>(p 2 1 k 2 )} never vanishes so that the functional traces on the RHS of 
( p.4|) are not suffering from any IR problem. In fact, if we write down the perturbative 
expansion of the traces, for instance, we see that all propagators are correctly cut off in 
the IR, and that loop momenta smaller than k are suppressed correctly. This would not 
have been the case if we had insisted on a positive Zk, setting Zk = —z k > 0. In this case 
the modified inverse propagator — |z/c|[p 2 — k 2 R^\p 2 / k 2 )], because of the relative minus 
sign between p 2 and the TZk-term, fails to suppress the IR modes. Even worse, it can 
introduce a spurious singularity at the value of p 2 for which p 2 — k 2 R^°\p 2 / k 2 ) = 0. 

At first sight the choice Z k = — z k > might have appeared to be the more natural one 
because only if Z k > the factor exp(— A k S) oc exp(— J 1Z k f 2 ) is a damped exponential 
which suppresses the low-momentum modes under the path integral. Nevertheless it was 
argued in QXJ] that the u Z k = 2fc-rule" is the correct choice both for z k > and z k < 0. The 
calculations in [ffj and all subsequent papers [EL [71 [53| were based upon this rule. On 
the one hand, this rule guarantees that the RG equation is well-defined and consistent. 
On the other hand, it is difficult to give a meaning to the Euclidean functional integral 
from which this RG equation was derived at the formal level. In the case Z k = z k < 
the factor exp(+ J \lZk\ V^ 2 ) is a growing exponential which seems to enhance rather than 
suppress the low-momentum modes. However, as suggested by the perturbative argument 
above, this conclusion is too naive probably. 

Let us now come to the case where the functions are not linear in D 2 but, say, 
of the form z { k 1] (D 2 ) 2 + z ( k ] D 2 . In this case the rule ( p.5|) demands that we choose the 
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corresponding operators (7Zk)ij in such a way that they contain cutoff terms adjusted to 
both the (D 2 ) 2 - and the D 2 -terms of the kinetic operator. For the above example, which 
is relevant to the _R 2 -truncation, (TZk)ij assumes the general form 



K t (p 2 ) = Zf 



2? fl<»> (£) + k* (> + * 2 A' ' i$) (3-7) 



where we omitted the indices referring to the types of fields. Obviously we have to 
set Zf = z« and zf = zjf 1 in order to achieve that the propagator and the cutoff 
combine as prescribed by the rule fl3.5|) . This leads to the modified inverse propagator 



z 



' ^\p 2 + k 2 R(°\p 2 /k 2 )] 2 + z^\p 2 + k 2 R ( - ">(p 2 /k 2 )]. For brevity we refer to this prescription, 



too, as the Z k = z k -ru\e. 

We believe that the Z\. = z^-rule is correct also for z k < 0, and that it is the relation 
between the manifestly well-defined flow equation and the formal path integral that needs 
to be understood better. Various attitudes are possible here. For instance, one could 
postulate that the fundamental definition of the theory is in terms of the flow equation 
rather than the path integral. Since the former is much better defined than the latter (in 
particular also with respect to the usual UV and IR problems) one would simply discard 
the path integral then. Another way out is to adopt the usual, albeit rather ad hoc, 
prescription of Wick-rotating the conformal factor (tp — * i(p) which turns the growing 
exponential into a decaying one. 

A much more attractive and less radical possibility is the following. Presumably it 
will be possible to construct an effective average action for Lorentzian quantum gravity by 
invoking a kind of stationary phase argument for the mode suppression. Then one deals 
with oscillating exponentials exp(iS) exp(iAfcS'), and apart from the trivial substitutions 
Tk — > — ir^, TZk — > — iR*k the flow equation remains the same as in Euclidean gravity. 
For Zk = Zk it has all the desired features, and Zk < poses no special problem for 
the path integral. It is interesting that there are also recent indications |38|] coming 



from the dynamical triangulation approach to quantum gravity which suggest that 



the Lorentzian path integral might have a better chance of being well-defined than the 
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Euclidean one. 

As the last possibility we mention the best of all situations, namely that in an exact 
treatment there are simply no factors Zk < 0. If this is actually the case, the conformal 
factor problem which we encounter in the Einstein-Hilbert and similar truncations would 
have the status of an unphysical truncation artifact. If so, the u Zk = Zfc-rule" could be 
interpreted as a device which helps in approximating as well as possible the exact RG 
flow by a truncated flow. 

It is one of the main results of the present paper that this scenario is indeed realized 
to some extent. We shall see that within the i? 2 -truncation those terms of which 
dominate at sufficiently large momenta have Zk > at least for large enough values of k 
(k ^> mpi). For too low scales (k < mpi) some of the z^s might turn negative, but at 
these scales the i? 2 -truncation becomes unreliable probably || so that the negative z^'s 
might be due to an insufficient truncation. It is not excluded that in the exact theory the 
dominating z^s are positive down to k = 0.Q 

As all those ZkS which determine the sign of the dominating contributions to are 
positive for large values of k, the cutoff has the standard suppression properties and 
is not plagued by any conformal factor problem for k — ■> oo. This, then, allows for an 
unambiguous investigation of the UV fixed point and its properties. It is quite remarkable 
that, as we shall see, all results concerning the fixed point are basically the same for the 
i? 2 -truncation (where Zk — z^ > at least for the dominating terms) and the Einstein- 
Hilbert truncation with both positive and negative factors Zf. = Zk- This is certainly a 
quite impressive confirmation of the Zk = z^-rule. 

3.3 The cutoff adapted to the i? 2 -truncation 

In the next section we shall see in detail that for the truncation studied in this paper 
we can comply with the Zk = Zk-rule by using the following cutoff operators for the 

4 Recent investigations in a scalar toy model jl0| indeed suggest that the conformal factor problem 
could be solved dynamically by strong instability-driven renormalization effects. 
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component fields: 
IT? 

y'^k) h T h T 



L>L> 



i + r ? r <*) { z f *V + yf~ hT R, k 2 R^(-D 2 /k 2 )} , 

Zfn 2 r v k 2 R {Q) {-D 2 /k 2 ) , 

X™ {-ID 2 k 2 R^\-D 2 /k 2 ) + k 4 R(°\-D 2 /k 2 ) 2 ) 

+ l -{z-fK 2 + y-fR,k 2 R^\-D 2 /k 2 )} , 



ft 



P* + j^jD^D^-D?)- 1 ) P k 



+D 2 J{D^ + ^- i D^D h 



+ (yfR + zf K 2 



Pk + -^-jD^Dvi-Dzyi ) P k 



Here Pi. is defined as 



Xf {-2D 2 k 2 R^ {-D 2 /k 2 ) + A; 4 P(°) (-£> 2 /k 2 ) 2 ) 



+ ^ |zf « 2 + yf R, k 2 R^\-D 2 /k 2 ) 



- (K k )^ T = zf vT r v k 2 R {0) {-D 2 /k 2 ) 
-{n k ) e5 = Zfk 2 R^{-D 2 /k 2 ) . 

P fe = _ J D2 + A ,2 jR (0) ( _ jD 2 /A; 2 ) 



(3.f 



(3.9) 



and the curly brackets denote the anticommutator, i.e. {A, B} = AB + BA for arbitrary 
operators A, B. The remaining cutoff operators which appear in Q2.9Q but are not listed 
in eq. (|3.8|) are set to zero. 

The constants X k , y k , and Z k will be adjusted later. It should be noted that the terms 
proportional to the X k s and 34's provide the cutoff for those contributions to 1^ which 



22 



come from the higher-derivative terms. For y]f^ T = y% a = y^ a = = X£ a = = 
= 0, eq. ( |3.8| ) actually boils down to the cutoff of type B used in in the context of 
the Einstein-Hilbert truncation. (This cutoff type has to be distinguished from the cutoff 
of type A used in the original paper fl[], which is formulated in terms of the complete 
fields and does not involve the component fields.) 

Each cutoff contains some "shape function" R(°\ A particularly suitable choice is the 
exponential shape function 

R^(y)=y[eMy)-ir 1 • (3-10) 

In order to check the scheme independence of universal quantities we employ a one- 
parameter generalization of ( p.lOj ), the class of exponential shape functions, 

R ( - \y;s) = sy[exp(sy)-ir 1 , (3.11) 

with the "shape parameter" s > parametrizing the profile of [35]. Another admis- 



sible choice we are going to use is the following class of shape functions with compact 
support: 



R {0 \y,b) = { 



1 y<b 

exp [(y- 1.5) _1 exp [(b - y)' 1 ]] b < y < 1.5 • (3.12) 
2/ > 1.5 



Here it is the shape parameter b G [0, 1.5) which parametrizes the profile of R^ [EI]. 



4 The i? 2 -truncation 



4.1 The ansatz 



In all previous papers [Q, |2(| |35|, |4], ^[ ||] the flow equation of QEG was used in the 
Einstein-Hilbert truncation. In this section we generalize this truncation by taking also 
an i? 2 -term with associated running coupling /3& into account and we derive the RG flow 
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within this "i? 2 -truncation" . We assume that, at the UV scale k — > oo, gravity in d 
dimensions is described by the action 

h\9] = S\g) = jd d x^! [T ^(-R(g) + 2X)+PR\g^ . (4.1) 

It consists of the conventional Einstein-Hilbert action and a higher-derivative term with 
bare coupling /3. In order to study the RG flow of T k [g,g] towards smaller scales k < k 
we employ a truncated action functional of the following form: 

r k [g,g] = J d d x^{2n 2 Z Nk (-R(g) + 2X k )+p k R 2 (g)} 

W^f J d d x sTgT v O^V) • (4.2) 

The ansatz (|4.2|) is obtained from S + by replacing 

G^G k = Z-\G, X^X k , ,9 -ft, a^Z N \a (4.3) 

so that its form agrees with that of the gravitational sector in the ansatz ( |3.1| ) with 

f k [g, g] = k 2 ^^± J d d x yfigT {H°g P *) • (4-4) 

In principle, also the gauge fixing parameter a should be treated as a scale-dependent 
quantity: a — > a k . Its evolution is neglected here for simplicity. However, setting a = 
by hand mimics a dynamical treatment of the gauge fixing parameter since a = can be 



argued to be a RG fixed point [42 



4.2 Projecting the flow equation 

The ansatz (|4.2p comprises three /c-dependent couplings. They satisfy the initial condi- 
tions A^ = A, Z N t = 1 which implies Gt — G, and (3t — (3. Here the UV scale k is taken 
to be large but finite. In order to determine the evolution of X k , Z^ k and (5 k towards 
smaller scales we have to project the flow equation onto the space spanned by the oper- 
ators J d d x J~g, J d d x JljR and J d d x ^R 2 ■ After having inserted the ansatz (|4.2|) into 
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(2) 

both sides of the flow equation and having performed the ^^-derivatives implicit in T k 
we may set = g^ u . As a consequence, the gauge fixing term drops out from the LHS 
which then reads 

d t T k [g,g] = J d d x^g{2n 2 [-R{g) d t Z Nk + 2d t (Z Nk \ k )] + R 2 (g) d t p k ) . (4.5) 

Obviously the LHS is spanned by the operators J d d x ^/g, j d d x ^R and J d d x ^fgR 2 ■ 
This means that we have to perform a derivative expansion on the RHS in order to extract 
precisely those contributions from the traces which are proportional to these operators. 
By equating the result to ( |4.5|) and comparing the coefficients we can read off the system 
of coupled differential equations for X k , Z^ k and f3 k . 

In order to make these technically rather involved calculations feasible we may insert 
any metric g^ v that is general enough to admit a unique identification of the operators 
spanning the truncated theory space. We exploit this freedom by assuming that g^ v 
corresponds to a maximally symmetric space. Such spaces form a special class of Einstein 
spaces and are characterized by 

- R _ _ . — R_ 

-tifivpa = ~^J^ gfia-gup) i tipu = ~~^9l lu \^'®J 

with the curvature scalar R considered a constant number rather than a functional of 
the metric. It is sufficient to employ spaces with R > 0, i.e. <i-spheres S . They are 
parametrized by their radius r which is related to the curvature scalar and the volume in 
the usual way, 

R = ^^> jd d x^g = y^(^)i. (4.7) 

We emphasize that the beta-functions of X k , Zjy k , and j3 k do not depend on this choice 
for g^ v \ it is simply a technical trick without any physical meaning. In principle the 
beta- functions could be computed without any specification of g^. 

While this projection technique is capable of distinguishing J d d x y/g oc r d from both 
J d d x y^-R oc r d ~ 2 and j d d x ^R 2 oc r d ~ A , it cannot disentangle the three (curvature) 2 - 
invariants f d d x^/g~R 2 ) J d d x y/g~R 2 u , and f d d x ^fg~R 2 ^ upa which are all proportional to 
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r d ~ 4 . If one wants to project them out individually one has to insert non-maximally 
symmetric spaces, but then the evaluation of the functional traces on the RHS of (|3.4[ ) 
is a rather formidable problem with the present technology. In fact, this is one of the 
reasons for omitting the other two (curvature) 2 -invariants from our truncation ansatz. 

In ref. we discussed already the expansion of fields defined on spherical back- 
grounds. Both the classical and the quantum TT-component fields can be expanded in 
terms of transverse-traceless tensor harmonics T l ™, transverse vector harmonics T^ m , and 
scalar harmonics T lm . They form complete sets of orthogonal eigenfunctions with re- 
spect to the corresponding covariant Laplacians. We summarize the main results of |2j in 
appendix |CJ In particular, the expansions of h? u , <fi, C M , and their classical counter- 
parts can be read off from eq. ( |C3| ) , while the remaining component fields are expanded 
according to 

oo oo £>j(d,0) 

u*) = E E UTjrix), ^) = E E °i™T lm (x), 

1=2 m=l 1=2 m=l 

oo Diidfl) oo D^dfl) 

= E E VimT lm (x), f}(x)=Yl E Vi m T lm (x). (4.8) 

1=1 m=l 1=1 m=l 

Similar expansions hold for the associated classical fields (expectation values). 

Contrary to eq. ( |C3|) , the summations in eq. (|4.8| ) do not start at I — 1 for vectors and 
at I — for scalars, but at / = 2 for £ M and a, and at I = 1 for the scalar ghost fields. The 
modes omitted here are the Killing vectors (T^ =1 ' m ), the solutions of the scalar equation 
( p.6|) (T l=1,m ), and the constants (T l=0 ' m=1 ) . As we mentioned in section [2|, these modes 
do not correspond to fluctuations of or the Faddeev-Popov ghosts. 

As in ||, we decompose the quantum field <fi into a part <pi spanned by the same set 
of eigenfunctions as a, and a part O containing the contributions from the remaining 
modes: 

1 AWO) oo 

<j>(x) = M*)+M*), <M^) = E E <j>i™T lm (x), 0i(^)=E E <PimT lm (x)(4.9) 

1=0 m=l 1=2 m=l 
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Due to the orthogonality of the spherical harmonics, 0o is orthogonal to 0i and a: (0i, 0o) 
= (a, (j) ) = 0. This implies (0,0) = (0 , 0o) + (0i,0i) and (er, 0) = (cx, 0i). As a 
consequence, splitting according to eq. (|4.9| ) ensures that any nonzero bilinear cross 
term of the scalar fields is such that the scalars involved can be expanded in the same set 
of eigenfunctions. Of course, the same holds for the corresponding classical fields 0o and 



4.3 Inserting the ansatz into the RHS of the RG equation 

Let us now start with the evaluation of the RHS of eq. ( |3.4|) . After having inserted the 
truncation ansatz (|4.2|) we identify the two metrics g pv and g^. Therefore it is sufficient 



to calculate the operators (T k 2 \g, g] + 7Z 



and (S^[g,g\ +K k [g\) 1 at g pv = g pv . To 



this end we first expand the ansatz ( |4.2|) according to 

T k \g + h,g] = T k [g,g] + 0{h) + Tf^[h]g] + 0(¥) 
and concentrate on the part quadratic in h pu , i.e. r^ uad [/i; g~\ . This leads to 



vr A [h-g] 



d d x^fgh pv \ K 2 Z Nk 



(4.10) 



(4.11) 



1 

4 

1 — a 



(r v D p D a - 8£D v D p ) 



a 



\ {^r u 9 P . - R 2 + 2g^R (-R p(T + DpDp - g pa D 2 ) 

+2R (5»R>* p - R% a v - Z5 v a b»D p + 2g pa D»D u + 25 p l 5»D 2 - 5 v a D p D») + 26 v a D p RD" 
-36 p l 5»D x RD x + A5 u a D^RD p - Ag pa D^RD u + g^g pa D x RD x + 2W v R p 



•-pa 



-ART [D p D a - g pa D 2 ) + 2D»D U D P D U - 4g^D 2 D p D a + 2g^g pa (D 2 ) 2 



At this stage g^ is still arbitrary. In order to (partially) diagonalize this quadratic form 
we insert the family of S d background metrics into eq. ( |4.11| ) and decompose h pv according 
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to eq. (|2.3|). Then we apply eq. ( fl9|) to the classical field <fi to decompose it as well. This 
yields 

1 IV 



+p k (RD 2 + G T (d) R 2 ) 



Z Nk K 2 (-D 2 + A T (d) R - 2X k ) 
- h T»u 



-Z Nk K 2 (-D 2 + A v (d, a)R- 2a\ k ) + G v (d) (3 k R 2 



a 



C S 2(d, a) Z Nk n 2 (-D 2 + A S2 (d, a)R + B S2 (d, a) X k ) 



+P k (H s (d) (D 2 ) 2 - G S i(d) {2RD 2 + R 2 ) 



a 



+20i 



C S 2(d, a) Css(d, a) Z Nk n 2 + fa (-H s (d) D 2 + 2G sl (d) R) 



x\/-D 2 \ -D 2 - 



R 



a 



<^>e{</>o,0i} 



d-l 

C S2 (d, a) C sl (d, a) Z Nk n 2 (-D 2 + A si (d, a)R + B S i{d, a) X k ) 



+fa (H s (d) (D 2 ) 2 - G S 2(d)RD 2 + G S3 (d) R 



(4.12) 



Here the various A's, B's, C's and G's and Hs are functions of the dimensionality d and 
the gauge fixing parameter a. The explicit expressions for these coefficients are given in 
appendix p.l| . 

This partial diagonalization is performed in order to simplify the inversion of the op- 
erator T^\g, g\ + 7Z k [g]. In fact, this is the main reason for using the TT-decomposition 
( p.3|) and specifying a concrete background. Note that in the pure Einstein-Hilbert trun- 
cation it is only the term in ( [4.1 1|) which is proportional to 1 — a that gives rise to mixings 
between the traceless part of h^ v and <f) and therefore necessitates the complete decom- 
position (|2.3|) . For a = 1, a complete diagonalization can be achieved by merely splitting 
off the trace part [p], |2|. This has to be contrasted with the i? 2 -truncation where the 
higher-derivative term introduces additional mixings between the traceless part and 0. 
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These cross terms do not vanish for a — 1. Hence the complete decomposition of is 
necessary for a partial diagonalization even in the case a = 1. 

At the component field level the cross terms boil down to a purely scalar cr-0 mixing 
term that vanishes for the spherical harmonics T l =°> m=1 anc l j , J=i,m_ gi nce these modes 
contribute to <f), but not to er, we cannot directly invert the associated matrix differential 
operator ( (1^ )$,• ) _ _ . As a way out, we split according to eq. (|4.9|) into O 

\ / i,j£{h T ,£,a,<f>} 

and 0i. This has the effect that only mixings between the scalars a and 0i survive, which 
can be expanded in the same set of eigenfunctions T lm starting at I = 2. Hence the 
resulting matrix differential operator ((rj^)y) _ is invertible. However, it 

should be noted that this additional split of leads to a slightly modified flow equation 
since it affects the matrix structure of this operator. In fact, the summation in the 
gravitational sector of eq. Q3.4Q now runs over the set of fields {h T , £, O > cf, 4>i}, with 

(R-k)fo$ = {TlkOfafa = (R-k)w and (R>k)afa = 

As a next step we calculate the contributions from the ghost fields appearing on the 



RHS of eq. (|3.4j) . For this purpose we insert the family of spherical background spaces 
S d into S'gh and set g^ v = g^ v . Then we use eq. (|2.7|) to decompose the ghost fields. This 
yields 

,R~ 



S gh [0, v,v;g] = V2 J d d x 



-T 



-D 



R 
1 



—D 



d 



Q 



(4.13) 



From now on the bars are omitted from the metric, the curvature and the operators 
D 2 and P h . 

Before we can continue with our evaluation we have to specify the precise form of the 



cutoff operators. Adapting them to 1^ and of eqs. ( |4.12| ), (|4.13|) by applying the 
rule ( P75| ) leads precisely to the structure fl$.8|) with the following choices for the X k 's, 
y^s and Z^s: 



,(2) 



\>(7<J 



X^O = ^ = Hs{d) fa ^ yWW = _^ 



y™ = 2G sl (d)Pk, y, 



4>o<t>o 
k 



G S 2(d)P k , Z, 



WW 
k 



'Nk 
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2f = -Z Nk , Zf- ~° = C S2 (d } a)C S3 (d, a)Z Nk , Zf = C S 2(d, a)Z Nk 
a 

Zjtih = zjkh = CM a ) Csi (d, a)Z Nk , zf" T = Zf = V2 . 



(4.14) 



- - (2) 

Thus, for g^ u = g^, the nonvanishing entries of the matrix differential operators T k +1Z 



and Sf? + 1Z k take the form 



(Tf\g,g} +TZ, 



h T h T 



rf '[g, g } + n k 



r k 2) [g,g} + n h 



T k 2) [g,g}+n k 



r k 2) [ g , g ] + n h 



S ( £[g,g] + n f 



Z Nk K 2 (P k + A T (d) R - 2\ k ) + fa (-RP k + G T (d) R 2 ) 

-Z Nk n 2 (P k + A v (d, a)R- 2a\ k ) + G v (d) fa R 2 , 

Cs 2 (d, a) Z Nk K 2 (P k + A S2 (d, a)R + B S2 (d, a) \ k ) 
+fa(H s (d)P 2 + G S i(d) (2RP k -R 2 )) , 
T k 2) [g,g}+TZ k ) 

C S 2(d, a) C S3 (d, a) Z Nk n 2 + fa (H s (d) P k + 2G S i R) 



xVPk\ Pk~ 



R 



d-1 



T {2) [g,g]+n h 



C S 2(d, a) C S i(d, a) Z Nk n 2 (P k + A si (d, a)R + B si (d, a) \ k ) 
+fa(H s (d)P 2 + G S 2(d)RP k + Gs 3 (d)R 2 ) , 



S^[g,g}+n, 



V2 



Pk 



R 



(S®[g, g ]+K k ) = ~(s^[g,g]+n k ) = V2 

\ / QQ \ J QQ 



For notational simplicity we set ( S^?[0, v,v; g}) = ( S^[g, g] ) with ipi, ip 2 £ h- 

Now we are in a position to write down the RHS of the flow equation with g^ = g^ u . 
We shall denote it S k (R) in the following. Obviously we need the inverse operators 



P k -2 



7( 2 ) r 



R 
~d 



(4.15) 



(rjj. + 1Z k ) 1 and (S^ + lZ k ) 1 . This inversion is carried out in appendix [A.l| . Inserting 
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the inverse operators into Sk(R) leads to the somewhat complicated result 



XI 



+Tr 



(IT) 



S k (R) = Tr (2ST2) [ (P k + A T {d) R-2X k + {Z Nk K 2 )~ l p k (-RP k + G T {d) R 2 )) 
(U-{Z Nk K 2 )-^ k RT^ 

(P k + A v (d, a)R- 2aX k + G v (d) {Z Nk K 2 )-^ k R 2 ) _1 jV 

{ (c S2 (d, a) (P k + A S2 (d, a)R + B S2 (d, a) X k ) 

+ (Z Nk K 2 )~ 1 p k (H s (d) P 2 + G S i(d) (2RP k -R 2 ))) 
x (c S2 (d, a) C S i(d, a) (P k + A si (d, a)R + B S i(d, a) X k ) 
+ {Z Nk K 2 y 1 p k (H s (d) P 2 k + G S2 {d)RP k + GsM R 2 ) ) 

- (C s2 (d, a) C S3 (d, a) + {Z Nk K 2 y^ k (H s (d) P k + 2G SI R) ) * P k (p k - 

x{ (C s2 (d, a) (P k + A S2 (d, a)R + B S2 (d, a) X k ) 
+(Z Nk K 2 y 1 p k (H s (d) P 2 + G S i(d) (2RP k - R 2 )) ) 
x (c S2 (d, a) C sl (d, a)M+ {Z Nk K 2 y l p k (H s (d) T 2 + G S2 {d) R%) ) 
+ {c S2 (d, a) C S i(d, a) (P k + A si (d, a)R + B S i(d, a) X k ) 
+(Z Nk K 2 )~ 1 p k (H s (d) P 2 + G S2 (d) RP k + G S3 (d) R 2 ) ) 
'c S2 (d, a) M + {Z Nk K 2 Y^ k (H s (d) T 2 + 2G sl (d) R%) ) 



R 



d- 



+2{c S2 (d, a) C S3 (d, a) + (Z Nk K 2 )-^ k (H s (d) P k + 2G S1 R) ) ^P k \jp^ 



1 



x 



2Z]y k K, 2 



d t [H s (d)p k 



d 



Pk W [Pk + j—^ D, BP» D„(-D*)-i\ P k 



+ D 2 J(D*) 2 <J 



d 
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+ [2G S i(d) p k R + C S2 (d, a) C S3 (d, a) Z Nk n' 



x 



-2Tr 



(IT) 



2Tr' 



(o) 



P fc -2^] A/q 



+ 



2Z Nk n 



Z=0 



0)fe 2 (d, «) C 5 i(d, a) f A,(d, 0) + ^R^^d, 0)/k 2 ) 



+A S i(d,a)R + B sl (d,a)X k 

+ {Z Nk K 2 )-^ k (H s {d) (A,(d,0) + fc a i2W(A|(d ) 0)/* a )) a 
+G S2 (rf) i? (A,(d, 0) + fc 2 i?(°)(Ai(d, 0)/A; 2 )) + G S3 (rf) i? 2 

+h(H s {d) (2A l (d,0)k 2 R^(A l (d,0)/k 2 ) + k A R ( -°\A l (d,0)/k 2 ) 2 ) 



+G S2 (d)Rk 2 R^(A l (d,0)/k 2 )^ 



The new quantities A/", A/o, 7i and 7^ introduced in eq. (|4.16|) are defined as 
AT = {2Z Nk )- 1 d t [Z Nk k 2 R^\-D 2 /k 2 )] 



1 - 2^( fc ) 



k 2 R^\-D 2 /k 2 ) + D 2 RW(-D 2 /k 2 ) , 

A/o = 2- 1 ^[A; 2 J R( )(- J D 2 /A ; 2 )] =fc 2 i?(°)(- J D 2 /A; 2 ) + J D 2 i?W'(-^ 2 A 2 
%. = (2^ k y 1 d t [p k k 2 R^(-D 2 /k 2 )] 

k 2 R(°\-D 2 /k 2 ) + D 2 RW(-D 2 /k 2 ) , 

T 2 ee (2^)- 1 94^4-2D 2 fc 2 i?(°)(-D 2 /A; 2 ) + A; 4 J R( )(- J D 2 A 2 ) 2 )] 



1 " ^Vp(k) 



(4.16) 



2P fc 



k 2 R^\-D 2 /k 2 ) + D 2 RW(-D 2 /k 2 
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--^(Jfe) {-D 2 k 2 R^(-D 2 /k 2 ) + k A R ( -°\-D 2 /k 2 ) 2 ) . (4.17) 

Here 

VN (k) = -d t \nZ Nk (4.18) 

and 

V(3 (k) = -d t lnP k (4.19) 

are the anomalous dimensions of the operators J d d x ^fgR and J d d x^R 2 , respectively. 
Furthermore, the prime at R^ denotes the derivative with respect to the argument. 



In eq. ( 4 . 1 6|) we refined our notation concerning the primes at the traces. From 



now on one prime indicates that the mode corresponding to the lowest eigenvalue has 
to be excluded, while two primes indicate the subtraction of the contributions from the 
lowest two eigenvalues. The subscripts at the traces describe on which kind of field 
the operators under the traces act. We use the subscripts (0), (IT) and (2ST 2 ) for 
spin-0 fields, transverse spin-1 fields, and symmetric transverse traceless spin-2 fields, 
respectively. 

4.4 The system of flow equations for A/j, gk and (3k 

Next we derive the flow equations for the couplings. In order to make the rather compli- 
cated calculations feasible we are forced to work from now on in the technically convenient 
gauge a = 1. Here we merely summarize the main steps, the details of the calculation 
can be found in appendix |A.2[ 

By expanding Sk{R) where R oc r~ 2 with respect to r and evaluating the traces 
by means of heat kernel techniques we extract those pieces from the RHS of the flow 
equation, eq. ( |4.16| ), which are proportional to the appropriate powers of the radius, i.e. 



r d oc f d d x yfg, r d ~ 2 oc f d d x ^/gR, and r d ~ 4 oc J d d x yfgR 2 '. Then we equate the result to 
the LHS, eq. fl4.5|) , and compare the coefficients of the various powers of r. This leads to 
a system of coupled differential equations for A^, Z Nk and /3k- 
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In order to present it in a transparent manner we introduce the dimensionless running 
cosmological constant 

X k = k~ 2 X k , (4.20) 

the dimensionless running Newton constant 

g k = k d - 2 G k = k d - 2 Z^ k G , (4.21) 

and the dimensionless running _R 2 -coupling 

P k = k 4 ~ d p k . (4.22) 

G k = G/Z Nk denotes the dimensionful running Newton constant. 

In terms of the couplings \ k , g k , and f3 k , our final result for the 3-dimensional flow 
equation reads 

d t \ k = (3 x (X k ,g k ,p k ;d) 

= Ai(A fe , g k , p k ; d) + r] N (k) A 2 (X k , g k , p k ] d) + r}p{k) A 3 (X k , g k , p k ] d) , (4.23) 

d t g k = (3 g (X k , g k , (3 k ; d) = [d - 2 + rj N (k)} g k , (4.24) 

d t l3 k = (3 p {X k , g k , p k ] d) = [4 - d - Vp (k)} f3 k . (4.25) 
The anomalous dimensions are explicitly given by 

rj N (k) = rj N (X k , g k , f3 k ; d) = (4.26) 

Bi{\ k) g k , p k ] d) [p k + C 3 (X k , g k , p k ] d)} - Ci(A fc , g k , p k ] d)B 3 (X k , g k , p k ] d) 
9k [1 - 9kB 2 {X k , g k , /3 k ; d)} [f3 k + C 3 (X k , g k , f3 k ; d)} + g k C 2 (X k , g k , p k ; d)B 3 (X k , g k , f3 k ; d) 

and 

Vp(k) = r]p{X k , g k , f3 k ; d) = (4.27) 

CijXk, g k , (3 k ] d) [1 - g k B 2 (X k , g k , p k ; d)} + gMM, 9k, Pk] d)C 2 (X k , g k , fa d) 
[1 - 9kB 2 (X k , g k , p k ] d)] \J3 k + C 3 (X k , g k , Pk] d)] + g k C 2 (X k , g k , Pk] d)B 3 (X k , g k , Pk] d) ' 



34 



The three /3-f unctions (3 X , (3 and (3p contain the quantities A4, B i} d, i = 1,2,3, which 
are extremely complicated functions of the couplings and the dimensionality d. The 
explicit expressions for these coefficient functions can be found in appendix [B|. They 
contain the new threshold functions \& and \I> of eqs. ( |A29 ) and ( |A30| ) which functionally 
depend on R(°\ They generalize the familiar threshold functions $ and $ which occur in 
the Einstein- Hilbert truncation. The three /3-functions are one of the main results of this 
paper. 



5 The fixed points 

5.1 Fixed points, critical exponents, and nonperturbative renor- 
malizability 

Because of its complexity it is clearly impossible to solve the system of flow equations for 
Afc, gk and fik, eqs. fl4.23|) , ( f4.24| ) and ( f4.25[ ), exactly. Even a numerical solution would 
be a formidable task. However, it is possible to gain important information about the 
general structure of the RG flow by looking at its fixed point structure. 

Given a set of /3-functions corresponding to an arbitrary set of dimensionless essential 
couplings gi(k), it is often possible to predict their scale dependence for very small and/or 
very large scales k by investigating their fixed points. They are those points in the space 
spanned by the gj where all /3-functions vanish. Fixed points are characterized by their 
stability properties. A given eigendirection of the linearized flow is said to be UV or IR 
attractive (or stable) if, for k — > 00 or k — > 0, respectively, the trajectories are attracted 
towards the fixed point along this particular direction. The UV (IR) critical hypersurface 
<Suv (Sir) in the space of all couplings is defined to consist of all trajectories that run into 
a given fixed point for k —>■ 00 (k — *■ 0). 

In quantum field theory, fixed points play an important role in the modern approach 
to renormalization theory |12| . At a UV fixed point the infinite cutoff limit can be taken 
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in a controlled way, the theory can be renormalized nonperturbatively there. As for 
gravity, Weinberg argued that a theory described by a RG trajectory lying on a 
finite- dimensional UV critical hypersurface of some fixed point is presumably free from 
unphysical singularities. It is predictive since it depends only on a finite number of 
free (essential) parameters. In Weinberg's words, such a theory is asymptotically safe. 
Asymptotic safety has to be regarded as a generalized, nonperturbative version of renor- 
malizability. It covers the class of perturbatively renormalizable theories, whose infinite 
cutoff limit is taken at the Gaussian fixed point g*j = 0, as well as those perturbatively 
nonrenormalizable theories which are described by a RG trajectory on a finite-dimensional 
UV critical hypersurface of a non-Gaussian fixed point g*j 7^ and are nonperturbatively 



renormalizable therefore 23 



Let us now consider the component form of the exact RG equation, i.e the system of 
differential equations 

kd kgi (k)=(3 i (g) (5.1) 

for a set of dimensionless essential couplings g(k) = (gi(fc), . . . ,g n (fc)). In an exact treat- 
ment the number n is infinite; in a specific truncation it might be finite. We assume that 
g* is a fixed point of eq. floTTl ), i.e. /3j(g*) = for all i = 1, . . . , n. We linearize the RG 
flow about g* which leads to 

n 

kd kgl (k) =Y, B v (S#) -g-j) (5-2) 
3=1 

where = djf3 { (g^) are the entries of the stability matrix B = (-By). Diagonalizing 
B according to S^BS* = — diag(6*i, . . . , 9 n ), S = (V 1 , . . . , V n ), where V 1 is the right- 
eigenvector of B with eigenvalue —9j we have 

n 

X>i^ = -OiV! , J=l,...,n. (5.3) 



0i 

(5.4) 



The general solution to eq. ( |5.2j ) may be written as 



&(*o=g*+i>^ (t) 
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Here 

n 

C/ = J>-%gi(*o) (5.5) 
i=i 

are arbitrary real parameters and is a reference scale. 

Obviously a fixed point is UV attractive for a given trajectory (i.e. attractive for 
k — > oo) only if all its Cj corresponding to negative 8j < are set to zero. Therefore 
the dimensionality Ayy = dim(iSuv) of the UV critical hypersurface belonging to this 
particular fixed point equals the number of positive ^/'s. Conversely, for a trajectory 
where all Cj corresponding to positive 9[ are set to zero, g* is an IR attractive fixed point 
(approached in the limit k — > 0). As a consequence, the IR critical hypersurface iS>ir of 
a fixed point has a dimensionality Air = dim(5iR) which equals the number of negative 
0/s. 

In a slight abuse of language we shall refer to the #r's as the critical exponents. 

Strictly speaking, the solution ( |5.4| ) and its above interpretation is valid only in such 
cases where all eigenvalues — Qj are real, which is not guaranteed since the matrix B is not 
symmetric in general. If complex eigenvalues occur one has to consider complex C/s and 
to take the real part of eq. ( |5.4| ), see below. Then the real parts of the critical exponents 
determine which directions in coupling constant space are attractive or repulsive. 

At this point it is necessary to discuss the impact a change of the cutoff scheme has 
on the scaling behavior. Since the path integral for depends on the cutoff scheme, 
i.e. on the A^S chosen, it is clear that generically the fc-dependent couplings and their 
fixed point values are scheme dependent. Hence a variation of the cutoff scheme, i.e. 
of TZk, induces a change in the corresponding B-matrix. So one might naively expect 
that also its eigenvalues, the critical exponents, are scheme dependent. In fact, this 
is not the case. According to the general theory of critical phenomena and a recent 



reanalysis in the framework of the exact RG equations |^3| any variation of the cutoff 
scheme can be generated by a specific coordinate transformation in the space of couplings, 
with the cutoff held fixed. Such transformations leave the eigenvalues of the B-matrix 
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invariant, so that the critical behavior near the corresponding fixed point is universal. 
The positions of fixed points are scheme dependent but their (non) existence and the 
qualitative structure of the RG flow are ^^-independent features. Quantities like the #/s 
which are ^-independent are called universal. Their residual scheme dependence present 
in an approximate treatment (truncation, etc.) can be used in order to judge the quality 
of the approximation. A truncation can be considered reliable only if it predicts the same 
fixed point structure for all admissible choices of lZ k . 

In the context of the i? 2 - truncation the space of couplings is parametrized by gi = A, 
g 2 = g and g3 = (3. The /3-functions occurring in the three flow equations 

dAk = /3 A (A fc , g k , f3 k ) , d t g k = (3 g (\ k , g k , p k ) , d t (3 k = Pp(\ k , g k , f3 k ) (5.6) 



are given in eqs. ( |4.23|) , ( f4.24|) and (|4.25|) , respectively 



In ref. the fixed point structure of the pure Einstein-Hilbert truncation was investi- 
gated. In this case the /3-functions were found to have both a trivial zero at A* = g* = 0, 
referred to as the Gaussian fixed point, and a non-Gaussian fixed point at A* ^ 0, 
g* ^ 0. As we will see in subsection |5^, the Gaussian fixed point is not present any more 
in the generalized truncation. This has to be contrasted with the non-Gaussian fixed 
point which is found with the i? 2 -truncation, too. In subsection |5.3| we study its cutoff 
dependence and the cutoff dependence of the associated critical exponents employing the 
above /3-functions with the families of shape functions (|3.11|) or (|3.12p inserted. 



5.2 The fate of the Gaussian fixed point 

In this subsection we study the fate of the Gaussian fixed point found in the context of the 
pure Einstein-Hilbert truncation. In we investigated the 2-dimensional RG flow near 
this fixed point (X*,g*) = (0,0) and discussed its stability properties. It is an important 
question how the situation changes by enlarging the parameter space. 

Quite remarkably, we find that in the 3-dimensional X-g-fi— space of the i? 2 -truncation 
there is no Gaussian fixed point, i.e. (\,g,/3) = (0,0,0) is not a simultaneous zero of 
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all three /3-functions. While (3 X and j3 g vanish at the origin (A, <?,/3) = (0,0,0), setting 
Afc = 9k = in (3p leads to 



(3JQ,Q,p k ;d)=f3J0,Q,Q;d)= ld V& . 



(5.7) 



The nonzero constant 7^ is given by 



ld = (47T)-5 {/l3l(rf) ^ /2 _ 2 (0) + M<0 ^/2-l(0) + M<*) ^ /2 (0)} 



(5. 



Here the /i,'s are defined as in subsection D.2. In d — 4 dimensions 
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419 
1080 



(4tt) 



(5.9) 



is a universal quantity since $o(0) = = 2*^2(0) = 1 independently of the cutoff, see 

appendix |]. 

Although there is no fixed point at the origin of the parameter space it is nevertheless 
very interesting to study the RG flow in the vicinity of (A, g, (3) = (0, 0, 0). For simplicity 
we restrict our considerations to the case d > 2. Expanding (/3 A , f3 g1 (3^) about the origin 
we obtain instead of ( [5~2"| ) the inhomogeneous system 

3 

k d k gi (k) = ld 5 ly3 + M H Si(*0 • ( 5 - 10 ) 

The linearized renormalization group flow is governed by the Jacobi-Matrix M = (Mij), 
= (9j/3j(0, 0, 0; d), which takes the form 



M 



(5.11) 



-2 v d d 
d-2 
y q d r d 4 - d J 

Its entries follow from the expanded /3-functions ( p3|) of appendix [E]. Here is d , s d and r d 
are (^-dependent parameters defined as 



v d = (rf-3)(4vr) 1 -$i /2 (0), 



(5.12) 
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ft = (4tt)-# {Mi)WO) + M^)W) +M<0*J/a(°)} » ( 5 - 13 ) 



r d = -(4vr) 1 - d [/i8 7 (d) + M<0 ^/ 2 (0)] 



X 



^(d) ^ /2 _ 2 (o) + ~h 35 (d) $ 2 d/2 _ x (o) + -^(d) $ d/2 (o) 



(5.14) 



d/2-2\ u J ' g"35l"; *,i/2-n u i T ^"361"; 
+ [/131(d) $i /2 _ 2 (0) + /l 32 (d) (0) + ^3 (d) $ d/2 (0)] 

M<0 $S /2 _ 2 (0) + h i0 (d) $^(0) + V(d) $ 2 /2 (0) + /142(d) $ 3 /2+1 (0) + I 



At this point it should be noted that the sub-matrix (My)jj 6 {i )2 } coincides precisely 
with the stability matrix of the Gaussian fixed point which was calculated in || in the 
Einstein-Hilbert truncation. 

Diagonalizing the matrix ( |5.11| ) yields the (obviously universal) eigenvalues $i = —2, 
$2 = d — 2 and $3 = 4 — d which are associated with the eigenvectors 

V l = (l,0,ft/(d-6)) T , V 2 = (z/ d ,l,(^ d + r d )/(2(d-3))) T , V 3 = (0, 0, 1) T (5.15) 



Eq. ( fy]D is valid only for d ^ 3,6. In d = 3 we obtain = (1,0,-^/3), I/ 2 = 
V^ 3 = (0,0,1), and in the 6-dimensional case the eigenvectors are V 1 = V 3 = (0,0,1), 
V 2 = (u 6 , 1, (^6^6 + r 6)/6). Thus in both cases the space spanned by the eigenvectors in 
only 2-dimensional, i.e. they do not form a complete system. For all values of d, including 
d = 3 and d = 6, the solutions for A& and g% obtained from the linearized system (|5.10|) 
assume the following form: 

2 



A* 



9k 



(Afc - Vd9k ) 



kr 



Vd 9k 



d-2 



9ko 1 To 



d-2 



(5.16) 



Since the expanded /3-function (3 of eq. (|E3|) does not depend on \ k and j3^ up 
to terms of third order in the couplings we can easily calculate also the next-to-leading 
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approximation for near the origin. In terms of the dimensionful quantity Gk this 
improved solution reads 



Gk — Gk a [l — u d Gh {k d 2 — k d 2 )] 



(5.17) 



with 



d-2 



Bi(0,0,0;d) = (4tt) 1 -5 {^3(d)$i /2 -i(0) + Mrf)^/2(0)} (5-18) 



a d-dependent parameter. It agrees with the u>d defined in || in the context of the pure 
Einstein-Hilbert truncation. For k <C \cudGk \~ 1 ^ d ~ 2 ^ and with the reference scale fco = 
(which is admissible only for specific initial conditions of the cosmological constant) eq. 
( CT ) yields 



G k = G [1 - u d G k d ~ 2 + O (G 2 k 2 ^)] 



(5.19) 



For the dimensionful cosmological constant we obtain from eq. ( 5.16Q 



A fe = A feo + u d G ko (k d - k d ) . 



(5.20) 



Eqs. ( |5.19j ) and ( |5.20| ) agree completely with the corresponding results from the Einstein- 
Hilbert truncation. 

Let us now discuss the solution for (3^- In order to derive it we start by picking i = 3 
in ( |5.10| ) and rewrite the corresponding equation as 



d k [k d ~ 4 p k ] = k d - 5 [ ld + s d \k + r d g k ) 



(5.21) 



Then we insert the solutions for gk and A& of eq. ( |5.16| ) into (|5.21 ). The resulting 
differential equation may easily be solved. For d ^ 3,4, 6, the solution reads 



n Id (Afcp ~ Vd9k Q ) ( ko\ VdSd + r d ( k_ 

Pk d-4 + d-6 \k) 2(d-3) 9k ° \k 



d-2 



+ 



flho 



Id _ (Afcp - Z^fffcoKri _ "dSd + T d 

d-4 d-6 2(d-3) 



9k 



d-4 



(5.22) 
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The solutions in d = 3, d — 4, and d = 6 can be obtained from ( |5.22|) by a careful 
evaluation of the limits d — > 3, d — ► 4, and d — > 6, respectively. In the most interesting 
case of d = 4 dimensions this leads to the following solution: 

(Afc — vigko) ^ 



0k 



Pk + 

(Afc 



+ r 4 419 
o 0*o + T7^7T( 47r J 



2 \ 

The parameters appearing in eq. fl5.23p are 
1 



I/ 4 ?4 + T 4 
H o # fe o 
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In 



(5.23) 



v A = 



Air 



$5(o) 



Q = (4 



7T 



559 71 



432 36 



*i(0) 
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24 



r 4 = (4 



7T 



i-3 



419 
1080 



| *i(0) + ^(0) 



299 



559 71 ~ 2/ , 347 -o, / 

+ — $?(o) + — $2(0) 



1728 288 



288 



(5.24) 



Employing the exponential shape function with s — 1 and inserting the corresponding 
values of the $^(0)- and $^ (O)-integrals given in appendix [F] we obtain 



1/4 = C(3)/(2tt) « 0.19 , ? 4 = (4tt) 



-2 



559 
432 



278, , , 347, , , 
— ln(2) - — ln(3) 



24 



0.20 



2817356 + 25(474 + 13vr 2 )(-559 + 4590 ln(2) - 2082 ln(3)) 



0.0051 . (5.25) 



49766400tt 3 

Let us now analyze the RG flow near the origin of the parameter space. Strictly 
speaking our analysis even extends to all points of the \-g-/3— space which satisfy |A|, \g\ <C 
1 and \(3 | <C l/|g|. This is because in any of the three /3-functions all terms of second and 
higher orders in /3j- appear as products g%(3™ with n > m. 

Since $1 = — 2 < and, for d > 2, $ 2 = d — 2 > 0, A& of eq. ( |5.16| ) starts growing as 
soon as k falls below k Q and the linearization breaks down, unless the couplings run along 
a trajectory which satisfies 



A. 



v<igk 



Afc — v d Gk k° 



(5.26) 
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for sufficiently small values of k 0, with Gk given by eq. ( |5.f 7| ). In this case both A 



and gk approach A = g = in the limit k — > as long as \(3}.\ -C l/|yjfc| is satisfied as 
well. Since 3>^/ 2 (0) depends on the shape function R(°\ is not a universal quantity. 
Therefore the slope of the distinguished trajectories characterized by (|5.26 ) is not fixed 
in a universal manner. 



Eq. (|5.26| ) is exactly the condition for the "separatrix" found in || in the context of the 



4-dimensional Einstein-Hilbert truncation. In the terminology of || , the separatrix is the 
"type Ha trajectory" that interpolates between the Gaussian and the non-Gaussian fixed 
point of the Einstein-Hilbert truncation, thereby separating the region of trajectories with 
Afc^o ~~ > ~ 00 (type la) from those which hit the boundary of parameter space (A = 1/2) 
for some finite value of k (type Ilia). In FIG. |l](a) we depict the separatrix, a type la and 
a type Ilia trajectory in the vicinity of (A, g) = (0,0). The plot should be thought of as 
a projection from \-g-(3— space onto its (3 = 0— plane. 

For the separatrix, the equation (|5.19 ) for the running of Gk is valid down to k — 
because gk stays near the origin. For d ^ 2 the parameters ^/2-i(^) anc ^ ^4/2$) 
appearing in uo^ are scheme dependent, and Ud is nonuniversal. In the most interesting 
case of d = 4, 

^4 = ^[l3$l(0) + 79^(0)] . (5.27) 

Since $J(0) and $1(0) are positive for any admissible shape function^ we can infer from 
eq. ([5.27D that is positive. Thus, if we define QEG with vanishing renormalized 
cosmological constant to be the theory described by a trajectory in \-g-f3— space, whose 
A- and ^-coordinates follow the separatrix in the limit k — > 0, eq. (|5.17|) implies that 
QEG is antiscreening in the IR, i.e. Gk decreases as k increases. 

Up to now we investigated the (3 = 0-projection of the flow on \-g-(3— space. Next we 

discuss the linearized RG flow of the /3-component. 

5 Using the exponential shape function with s — 1, for instance, we have $1(0) = 7r 2 /6, ^(O) = 1 
so that k>4 w 1.33. Furthermore, we have &l(0) = 2£(3) where £ denotes the Riemann zeta function, and 
thus v A w 0.19. 
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5.2.1 d = 4 

In d = 4 the solution for j3 k , eq. (|5.23|) , diverges in the limit k — > 0. If 7^ fd<7fc as 
/c — >• 0, the leading divergence is quadratic in k. However, in this case the linearization 
cannot be trusted down to arbitrarily small values of k anyhow since trajectories with 
X k 7^ v d g k ultimately run away from (A, g) = (0, 0) for k — > 0. 

For the distinguished trajectories which satisfy the separatrix condition ( |5.26| ) for 
k — > 0, the coefficient of the term oc (ko/k) 2 in ( |5.23| ) vanishes and only a logarithmic 
running with a universal coefficient remains: 

V&A + 7"4 , 419 , t / k \ . 2\ 



A = ft, - 2 g fc0 + Y^( 47r ) ln ( H + °(* ) • ( 5 - 28 ) 

Since higher orders in f3 k appear exclusively as products g^fi™ with n > m, the vanish- 
ing of g k (ln(k/ko)) m oc (k / ko) 2n (ln(k / ko)) m in the limit k — >• then implies that terms of 
order /3| remain negligible as — > 0. As a consequence, the linearization does not break 
down for k — > although (3k diverges in this limit. 

According to eqs. ( |5.16| ), (|5.26|) , A^ and g k quickly approach A = g = so that the 



corresponding trajectories run almost along the /3-axis for k — > 0, and the RG flow becomes 
essentially one-dimensional. This logarithmic running of j3 k was expected on the basis of 
conventional perturbation theory ||31|| . We observe that \f3 k \ decreases logarithmically 
with increasing k. This is what is usually referred to as the "asymptotic freedom" of 
the (curvature) 2 -coupling. We emphasize, however, that according to our results this 
logarithmic running occurs only close to g = A = and does not represent the true short- 
distance behavior of the theory. In fact, we shall find that (3 k runs towards a fixed point 
value /3* for k — > 00. 

FIG. 0(b) shows three typical trajectories of the i? 2 -truncation close to (\,g,/3) = 
(0,0,0); all of them satisfy the separatrix condition (|5.26|) . Their /3-component diverges 
logarithmically towards — 00 as k goes to zero, which is due to the positive coefficient in 
front of the ln(/c/fco)-term in ( |5.28| ). In this figure we also depict the common projection 
of the trajectories onto the X-g— plane. It coincides precisely with the separatrix of the 
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(a) (b) 

Figure 1: The case d = 4: (a) Type la, type Ha, and type Ilia trajectories (from left 
to right) obtained from the X-g— projection Q5.16Q . They coincide precisely with the 
corresponding trajectories of the Einstein-Hilbert truncation. The type Ha trajectory is 
the separatrix with Xk=o = 0, which separates the region of trajectories with Afc_»o — > — oo 
(type la) from those running towards more positive A's (type Ilia), (b) Three typical 
trajectories of the linearized X-g- [3— equation. They correspond to different values of j3k , 
but all of them satisfy ( |5.26| ). We also depict their projection onto the X-g— plane. It 
coincides with the separatrix of the pure Einstein-Hilbert truncation. 

Einstein-Hilbert truncation i.e. the curve in FIG. [T](a) that hits (A, g) = (0,0). 
Conversely, all trajectories of the i? 2 -truncation satisfying (|5.26 ) represent specific "lifts" 



of the separatrix with nonvanishing f3- components; they are distinguished by their j3k - 
and <7fc - values. 

5.2.2 

Let us now discuss the /3-evolution for 2 < d ^ 4. Again, the linearization breaks down 
for trajectories which do not satisfy ( |5.26| ) for k — > since |Afc_»o| ^ oo in this case. 



Therefore we restrict our considerations to the trajectories with A& = v^gu for sufficiently 
small k. In this case the second, quadratically divergent term of eq. (|5.22|) drops out, and 
the only powers of k which occur in f3k are k d ~ 2 and k A ~ d . 

In 2 < d < 4, both $2 = d — 2 and ^3 = 4 — d are positive which implies that the 
RG trajectories considered are attracted towards (X,g,/3) = (0,0,7^/(^ — 4)) as k is sent 
to zero. 



45 



For d > 4, $3 is negative and thus /3& contains a divergent term oc k A ~ d . As a con- 
sequence, the coefficient of this term must vanish, if a trajectory is to hit the point 
(X,g,(3) = (0, 0, 7d/(d — 4)) in the limit — >■ 0. The distinguished trajectory which runs 
into this point as k — > satisfies, for sufficiently small values of k, 

a v <&d + T d v d q d + T d 

Pk ~ 2(d-3) ^~ 2(d-3)u d k 

^ p k = ^g±^ G k k 2 ^ = V f* + l d \ k k d ~ 6 . (5.29) 
H 2(d-3) 2(d-3)z/ d v 7 

For all other trajectories the /^-component diverges for k —>■ 0. However, higher orders 
of are again suppressed by powers of g k and may therefore be neglected. (Note that 
lim^^o dkPk 1 li m fc-*o /c n ( d_2 ) _m ( d ~ 4 ) = for n > m.) As a consequence, the linearization 
can be trusted down to arbitrarily small scales k even in this case. The shape of the 
corresponding trajectories resembles the one found in d = 4. While \/3k\ — > oo, the A- and 
(^-components approach A = g = in the limit k — > 0. Thus, for sufficiently small fc, the 
trajectories are almost straight lines which virtually coincide with the /3-axis. 

Having a closer look at the /3-functions one recognizes that the IR scaling behavior in 
d 7^ 4 dimensions is actually governed by a "quasi-Gaussian" fixed point at 

(K,g*,P*) = (0,0,7 d /(d-4)) . (5.30) 

The quasi-Gaussian fixed point is not present in d — 4. Linearizing the RG flow about 
this fixed point yields essentially the same results as our expansion about (0, 0, 0) above. 
The linearized /3-functions with stability matrix B, and the linear solutions associated to 
this fixed point, may be obtained from eqs. (|5.10| ), ( |5.11|) , (|5.16| ), (|5.22|) and (|5.29|) simply 



by replacing r d with 



2(4 



r d = r d - " V J^ 4 (/i 39 (d) + h i5 (d) ^ /2 (0)) . (5.31) 

In particular, B = lSA(r d — > r d ). The constants Qj = — i9j assume the meaning of critical 
exponents now, and their signs determine the dimensionality Air of the (truncated) IR 
critical hypersurface (Sir of the quasi-Gaussian fixed point. 



,1-d 
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In 2 < d < 4 we have one positive critical exponent Q\ > and two negative critical 
exponents 9 2 ,0 3 < 0. Therefore, within the truncation, A IR = 2, as suggested by the 
corresponding solutions discussed above. 

In d > 4, 6i and #3 are positive and 62 is negative. Hence, in this case Air = 1, i.e. «Sir 
consists of a single trajectory. For sufficiently small values of k this IR critical trajectory 
is given by eq. ( |5.29| ) with r d replaced with r d . Since the parameters u d , q d and r d contain 
.R^-dependent integrals $^(0), $^(0), they are not universal. Therefore the slopes in 
both directions of the distinguished trajectory ( |5.29| ) are not fixed in a universal manner. 
This is in accordance with the general expectation that the eigenvalues of B should be 
universal, but not its eigenvectors. 




(a) (b) 

Figure 2: (a) The case d = 3: three typical trajectories of the linearized flow. They 
correspond to different values of j3k , but all of them satisfy ( |5.26| ). As a consequence, all 
3 trajectories hit the quasi- Gaussian fixed point for k — > 0. (b) The case d = 5: three 
typical trajectories of the linearized flow corresponding to different values of fik - All of 
them satisfy ( |5.26| ), but in contrast to the other two curves, the one in the middle satisfies 
also (|5.29| ). As a consequence, it hits the quasi- Gaussian fixed point for k — > 0. In both 
(a) and (b) we also depict the projection of the curves onto the (3 = 0— plane. 

We illustrate our results for d 7^ 4 in FIG. 0. In FIG. 0(a) we consider d = 3 and 
in FIG. 0(b) the 5-dimensional case. Each figure shows three typical trajectories in the 
vicinity of the quasi- Gaussian fixed point. All of them satisfy ( |5.26|) , so that in both d = 3 
and d = 5 the projections of the 3 trajectories onto the (3 = 0— plane coincide with the 
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separatrix. 

In d = 3, FIG. 0(a), all 3 trajectories hit the fixed point, independently of their /3k - 
value. As we already pointed out above, the quasi- Gaussian fixed point in d < 4 is IR 
attractive for all trajectories satisfying (|5.26| ). In d > 4 this is no longer the case. This is 
confirmed by FIG. 0(b) for d — 5. Here only one of the trajectories hits the quasi- Gaussian 
fixed point for k — > 0, and this is precisely the one which satisfies the additional condition 
( |5.29| ). The other two trajectories shown in FIG. @(b) correspond to /^-values which are 



different from the one in ( |5.29| ) and thus their /^-component diverges in the limit k — > 0. 
Depending on the /3fc -value, /3k runs towards +oo or — oo. 

5.3 The non-Gaussian fixed point 

Now we turn to the nontrivial simultaneous zeros of the set of /3-functions {(3\, (3 g , /3g} 
given by eqs. ( |4.23[ ), Q4.24| ), ( |4.25[ ). Such non-Gaussian fixed points with A*, g*, (3* all 



different from zero have the anomalous dimensions 

= 2 — d , 77,3* = d — 4 (5.32) 
which follow immediately from eqs. Q4.24 ) and ( |4.25|) . 



5.3.1 Results obtained from the pure Einstein-Hilbert truncation 

In d = 4 dimensions, and for the cutoff of the type A introduced in || , the non-Gaussian 
fixed point of the pure Einstein-Hilbert truncation was first discussed in [20, ||, and in 
ref. |35| the a- and -dependence of its projection (0, g*) onto the (^-direction has been 
investigated. However, since for a ^ 1 the cutoff of type A was defined in [34| by an 



ad hoc modification of the standard one-loop determinants it is not clear whether it can 
be derived from an action A^S, except for the case a — 1 [1J. Since a specification of 
A^S is indispensable for the actual construction of T^, the status of the results derived in 
]35] is somewhat unclear. In ref. we performed a comprehensive analysis of the fixed 
point properties using different cutoffs of type B, for which a AkS is known to exist. In 
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particular, we investigated the cutoff scheme dependence of various universal quantities of 
interest, both by looking at their dependence on the shape function and by comparing 
the "type A" and "type B" results. 

In this respect universal quantities are of special importance because, by definition, 
they are strictly cutoff scheme independent in the exact theory. Any truncation leads to 
a scheme dependence of these quantities whose magnitude is a measure for the reliability 
of the truncation j|3| . Typical examples of universal quantities are the critical exponents 
6j. The existence or nonexistence of a fixed point is also a universal, scheme independent 
feature, but its precise location in parameter space is scheme dependent. Nevertheless 
it can be argued that, in d = 4, the product g*A* is universal JP1 [| while and A* 
separately are not. 

For later comparison with the i? 2 -truncation, let us briefly list some of the results we 
obtained in with the pure Einstein-Hilbert truncation: 

(Ie.h.) Universal Existence : Both for type A and type B cutoffs the non-Gaussian 
fixed point exists for all shape functions we considered. This result is highly nontrivial 
since in higher dimensions (d > 5) the fixed point exists for some but does not exist for 
other cutoffs ||. 

(2e.h.) Positive Newton Constant : While the position of the fixed point is scheme 
dependent, all cutoffs yield positive values of g* and A* . A negative g* might be problematic 
for stability reasons, but there is no mechanism in the flow equation which would exclude 
it on general grounds. 

(3e.h.) Stability: For any cutoff employed, the non-Gaussian fixed point is found to be 
UV attractive in both directions of the X-g— plane. Linearizing the flow equation according 



to eq. ( |5.2|) we obtain a pair of complex conjugate critical exponents 6\ = Q\ with positive 
real part 9' and imaginary parts ±8". Due to the positivity of 6', all trajectories in its 
basin of attraction hit the fixed point as k is sent to infinity. Because of the nonvanishing 
imaginary part 9" the trajectories spiral into the fixed point for k — > oo. 

Solving the full, nonlinear flow equations || shows that the asymptotic scaling region 
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where the linearization is valid extends from k u = "oo down to about k ~ mpi with the 

— 1/2 

Planck mass defined as mp\ = Go ' . It is the regime above the Planck scale where the 
asymptotic freedom of Gk sets in. 

(4e.h.) Scheme- and Gauge Dependence: The critical exponents are reasonably con- 
stant within about a factor of 2. For the gauges a = 1 and a = 0, for instance, they 
assume values in the ranges 1.4 < & < 1.8, 2.3 < 0" < 4 and 1.7 < < 2.1, 2.5 < 0" < 5, 
respectively. The universality properties of the product A* are much more impressive 
though. Despite the rather strong scheme dependence of g* and A* separately, their prod- 
uct exhibits almost no visible -R^-dependence. Its value is g*A* ~ 0.12 for a = 1 and 
g*\* m 0.14 for a = 0. The differences between the "physical" (fixed point) value of the 
gauge parameter, a = 0, and the technically more convenient a = 1 are at the level of 
about 10 to 20 per-cent. 

The above results suggest that the UV attractive non-Gaussian fixed point occuring 
in the Einstein-Hilbert truncation is very unlikely to be an artifact of this truncation but 
should rather be the projection of a fixed point in the exact theory. We interpreted them 
as nontrivial indications supporting the conjecture that 4-dimensional QEG is "asymp- 
totically safe" in Weinberg's sense. 

5.3.2 Results obtained from the i? 2 -truncation 

The actual justification of a truncation is that when one adds further terms to it its 
physical predictions do not change significantly any more. In order to test the stability 
of the Einstein-Hilbert truncation against the inclusion of other invariants we shall now 



reanalyze the non-Gaussian fixed point in the generalized truncation Q4.2D including the 
i? 2 -term. Starting from the (3- functions of the i? 2 -truncation, eqs. ( 4.23 ), (|4.24|) and 
( |4.25| ), we determine the location of the fixed point in \-g-(3— space and the linearized 
flow in its vicinity. Then we investigate the residual cutoff scheme dependence of the 
associated universal quantities, and we compare our results to those obtained from the 
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pure Einstein-Hilbert truncation. 

Note that, contrary to the pure Einstein-Hilbert truncation, only a cutoff of type B is 
used in the context of the generalized truncation. Therefore we omit the specification of 
the cutoff type when we refer to results obtained from the i? 2 -truncation. 

Location of the fixed point (d = 4) 

In a first attempt at finding the non- Gaussian fixed point in the i? 2 -truncation we neglect 
the cosmological constant and the coupling of the i? 2 -invariant. We approximate A& = 
A* = 0, (3k = (3* = 0, thereby projecting the renormalization group flow onto the one- 
dimensional space parametrized by g. In this case the non-Gaussian fixed point is obtained 
as the nontrivial solution of (3 g (0, g*, 0; d) = 0. It is determined in appendix [E| with the 
result given by eq. (|E2|). For any d, this solution coincides precisely with the analogous 
approximate solution (H2) of ref. || with a = 1, obtained in the pure Einstein-Hilbert 
truncation. In order to get a numerical value for the fixed point we have to specify R^°\ 
Inserting the exponential shape function with s = 1 into eq. ( E2 ) and setting d = 4 leads 
to g* » 0.590. 

Assuming that for the combined X-g-(3 system the numbers A*, g* and /3* are of the 
same order of magnitude as g* above we expand the /3-functions about Qki flk) — 
(0, 0, 0) and neglect terms of higher orders in the couplings. Again in appendix [E] we 
determine the non-Gaussian fixed point from the corresponding system of differential 
equations. Inserting the shape function ( |3.1(J| ) and setting d = 4, we find (\*,g*,(3*) ~ 
(0.287,0.751,0.002). Quite remarkably, for any cutoff A* and g* agree perfectly with 
the corresponding values obtained in by the same approximation applied to the pure 
Einstein-Hilbert truncation. 

In order to determine the exact position of the non-Gaussian fixed point (A*, g*, (3*) 
we have to resort to numerical methods. Given a starting value for the fixed point, 
for instance one of the approximate solutions above, the program we use determines a 
numerical solution which is exact up to an arbitrary degree of accuracy. Under the same 



51 



conditions as above, i.e. for s = 1 and d = 4, we obtain 



(A*, g„ A,) = (0.330, 0.292, 0.005) (5.33) 

In the pure Einstein-Hilbert truncation the corresponding coordinates of the fixed 
point are (A*,<7*) = (0.348,0.272) 0. Obviously the values of A* and g* are almost the 
same in both cases. While A* and g* are of the same order of magnitude, we find that /3* 
is significantly smaller than A* and g*. 

In order to test whether these properties of the fixed point coordinates are univer- 
sal we study their scheme dependence by looking at the s- or 6-dependence introduced 
via the one-parameter families of shape functions (|3.11| ) or (|3.12|) , respectively. Here s 



parametrizes the family of exponential shape functions ( |3.11| ), while the shape parameter 



b allows us to change the profile of the shape functions with compact support ( p3.12| ). 



As for the family of exponential shape functions, we are forced to restrict our consid- 
erations to shape parameters s > 1. This is because for s < 1 the numerical integrations 
are plagued by convergence problems. They are due to the fact that in d = 4 some of the 
threshold functions appearing in (3 X , f3 g and /3g diverge in the limit s — > 0, see also |3~5 |. 



As for the family of shape functions with compact support, we have to restrict ourselves 
to b < 1.2 for similar reasons. Here R^(y;b) approaches a sharp cutoff as b — > 1.5, 
which introduces discontinuities into the integrands of the threshold functions $^ and 
Already for b > 1.2 the /3-functions start to "feel" the sharp cutoff limit, which leads to 
convergence problems. 

As in the case of the pure Einstein-Hilbert truncation our results establish the 
existence of the non-Gaussian fixed point in a wide range of s- and b- values. As expected, 
the position of the fixed point turns out to depend on s or b, i.e. on the cutoff scheme, 
but the crucial point is that it exists for any of the cutoffs employed. FIGS. || and § show 
its coordinates (A*, g*, /?*) as well as the product g*A* for the shape functions ( |3.11| ) and 
( |3.12| ), respectively. In Fig|^(a) we plotted the various quantities in the range 1 < s < 5 
where the largest changes in A* and g* occur, but we calculated them for 1 < s < 30. 
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(a) (b) 

Figure 3: (a) g*, A*, and g^A* as functions of s for 1 < s < 5, and (b) /3* as a function of 
s for 1 < s < 30, using the family of exponential shape functions. 

For every shape parameter s or b, the values of A* and g* are almost the same as those 
obtained with the Einstein- Hilbert truncation [[|. As a consequence, the product g*A* is 
again almost constant and its value differs only slighlty from the one in || for the same 
gauge a = 1. Both FIGS. |](a) an d H(a) suggest the universal value g*A* ~ 0.14 while 
we obtained g^A* ~ 0.12 from the pure Einstein-Hilbert truncation. Thus we may expect 
that our g^-value is precise at the 10 to 20 per-cent level. Presumably this degree of 
precision is the best we can achieve in the present calculation because we saw already 
that the error due to using a = 1 instead of the "correct" a = leads to an uncertainty 
of the same size. 

Furthermore, our results show that /3* is always significantly smaller than g* and A* for 
both families of shape functions, which is quite remarkable. Within the limited precision 
of our calculation this means that in the three-dimensional \-g-f3— space the fixed point 
practically lies in the A-g— plane with (3 = 0, i.e. on the parameter space of the pure 
Einstein-Hilbert truncation. 

It is also interesting to note that the scheme dependence of /3* is unexpectedly small. 
As for the family of exponential shape functions ( |3.11| ), the function depicted in FIG. 
|3](b) develops a plateau-like shape for not too small values of s. Employing the family of 
shape functions with compact support, the scheme dependence of /3* is even weaker. The 
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Figure 4: (a) g*, A* and g*\*, and (b) /3* as functions of b for < b < 1.2, using the 
family of shape functions with compact support. 

function f3*{b) plotted in FIG. §(b) is almost constant in the range < b < 1.2. More- 
over, the positions of the two plateaus are nearly identical. While FIG. ||(b) suggests the 
value & « 0.0031, we obtain & « 0.0036 from FIG. |(b). This indicates that in d = 4 
dimensions also /3* might be a universal quantity. 

The linearized flow (d = 4) 

Let us now analyze the critical behavior near the non-Gaussian fixed point. Quite remark- 
ably, the non-Gaussian fixed point of the i? 2 -truncation proves to be UV attractive in any 
of the three directions of \-g-f3— space, for all cutoffs used. The linearized flow in its vicin- 
ity is always governed by a pair of complex conjugate critical exponents 9\ = 9' + \9" = 9\ 
with 9' > and a single real, positive critical exponent 63 > 0. (We define Q\ as the critical 
exponent with the positive imaginary part so that 9" > 0.) The general solution to the 
linearized flow equations is obtained by taking the real part of eq. (|5.4j) . Introducing the 
RG time t = \n(k/ko) it may be written as 

(h,g k ,Pk) T = (A*,3„&) T + 2j [ReC cos(0"t) +ImC sin Re V 

+ [ReC sin (0"t)-ImC cos (9" t)) Im V j e -0 '* + C 3 V 3 e - 93i (5.34) 
with arbitrary complex C = C\ = (C2)* an d arbitrary real C3. Furthermore, V = V 1 = 
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(V 2 )* and V 3 are the right-eigenvectors of the stability matrix (-3y)ij'e{A,5,#} with eigen- 
values —Q\ = —9\ and —#3, respectively. Obviously the conditions for UV stability are 
& > and 63 > 0. They are indeed satisfied for all cutoffs. As a consequence, all RG 
trajectories which reach its basin of attraction hit the fixed point as t is sent to infinity. 
The trajectories (|5.34Q comprise three independent normal modes with amplitudes pro- 
portional to ReC, ImC and C3, respectively. The first two are of the spiral type, the 
third one is a straight line. 

Let us illustrate these features by means of an example. For the exponential shape 
function ( p.ll|) with s = 1, for instance, we have (A*, g*, /3*) = (0.330,0.292,0.005). The 
corresponding stability matrix B takes the form 

' 8.83 2.61 401.75 ^ 
B = - 6.18 4.46 89.24 . (5.35) 
, 0.29 0.32 19.82 j 

It leads to the pair of complex critical exponents 9\ = 9\ with 9' = 2.15, 9" = 3.79, and 
to the real critical exponent #3 = 28.8. For the associated right-eigenvectors we find 

ReV = (-0.164, 0.753, -0.008) T , 
ImV = (0.64,0,-0.01) T , 

V s = -(0.92, 0.39, 0.04) T . (5.36) 

(The vectors are normalized such that ||V|| = ||V 3 || = 1.) In FIG. [5] we show a typical 
trajectory which has all three normal modes excited with equal strength (ReC = ImC = 
l/y/2, C 3 = 1). All its way down from k " = " 00 to about k = mpi it is confined to a very 
thin box surrounding the (3 = 0— plane, i.e the parameter space of the Einstein-Hilbert 
truncation. 

In fact, the linearized flow is characterized by the following quite remarkable properties, 
independently of the cutoff. They all indicate that, close to the non-Gaussian fixed point, 
the RG flow is rather well approximated by the pure Einstein-Hilbert truncation. 
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(a) (b) 

Figure 5: Trajectory of the linearized flow equation obtained from the i? 2 -truncation for 
1 < t — \n(k/k ) < oo. In (b) we depict the eigendirections and the "box" to which the 
trajectory is confined. 

(a) The /3-components of ReV and ImV are very tiny. Hence these two vectors 
span a plane which virtually coincides with the X-g— subspace at (5 — 0, i.e. with the 
parameter space of the Einstein-Hilbert truncation. As a consequence, the ReC- and 
ImC- normal modes are essentially the same trajectories as the "old" normal modes 
already found without the i? 2 -term. Also the corresponding 9'- and ^"-values coincide 
within the scheme dependence, see below. 

(b) For all cutoffs employed, the new eigenvalue 9s introduced by the i? 2 -term is 
significantly larger than 9', see below. When a trajectory approaches the fixed point from 
below (t — > oo), the "old" normal modes oc ReC, ImC are proportional to exp(—9't), but 
the new one is proportional to exp(— 9 3 t), so that it decays much more quickly. For every 
trajectory running into the fixed point, i.e. for every set of constants (ReC, ImC, C 3 ), 
we find therefore that, once t is sufficiently large, the trajectory lies entirely in the Re V- 
Im V— subspace, i.e. the f3 = 0— plane practically. 

Due to the large value of 6*3, the new scaling field is very "relevant". However, when 
we start at the fixed point (t u = "00) and lower t it is only at the low energy scale k ~ mpi 
(t ~ 0) that exp(— 9 3 t) reaches unity, and only then, i.e. far away from the fixed point, 
the new scaling field starts growing rapidly. 
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(c) Since the matrix B is not symmetric its eigenvectors have no reason to be orthogo- 
nal. In fact, we find that V 3 lies almost in the Re V^-Im V— plane. For the angles between 
the eigenvectors given above we obtain <(ReV^, ImV^) = 102.3°, <(KeV, V 3 ) = 100.7°, 
<(lmV,V 3 ) = 156.7°. Their sum is 359.7° which confirms that Rel/, ImV^ and V 3 are 
almost coplanar. This implies that when we lower t and move away from the fixed point 
so that the V 3 - scaling field starts growing, it is again predominantly the J d d x yfg- and 
/ d d x y^-R-invariants which get excited, but not J d x yfgR 2 in the first place. 

Summarizing the three points above we can say that very close to the fixed point the 
RG flow seems to be essentially two-dimensional, and that this two-dimensional flow is 
well approximated by the RG equations of the Einstein-Hilbert truncation. 

Scheme dependence of the critical exponents (d = 4) 

As we pointed out already the critical exponents are universal in an exact treatment, but 
in a truncated parameter space a scheme dependence is expected to occur as an artifact of 
the truncation. We may use it to judge the quality of our truncation. Also in this respect 
the i? 2 -truncation yields satisfactory results, which we display in FIGS. ^| and [7]. FIGS. 
^|(a) and 0(a) show the real and the imaginary part 9' and 9" of the complex conjugate 
pair 9i = 9\ while #3 is depicted in FIGS.|6|(b) and 0(b). The plots in FIG. [] are based 
on the family of exponential shape functions ( |3.11| ) and those in FIG. |7] are obtained by 
employing the family of shape functions with compact support ( |3.12|) . They display the 
s- and the 6-dependence of the critical exponents, respectively. 

As for the complex conjugate pair of critical exponents, the scheme dependence is of 
the same order of magnitude as in the case of the Einstein-Hilbert truncation 0. While 
the scheme dependence of 9" is weaker than that found in we see that it is slightly 
larger for 9'. For the exponential shape functions with 1 < s < 30, 9' and 9" assume values 
in the ranges 2.1 < 9'{s) < 3.4 and 3.1 < 9"(s) < 4.3, respectively. Employing the shape 
functions with compact support leads to a weaker dependence on the shape parameter 
b. However, the corresponding values 9'(b) and 9"{b) are in good agreement with those 
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Figure 6: (a) 9' = Re#i and 9" = Im^, and (b) #3 as functions of s, using the family of 
exponential shape functions. 

obtained with the exponential cutoffs. In fact, they all lie in the 9'(s)- and 9"{s)— intervals 
given above. The average values of 9' and 9" are slightly larger than those obtained from 
the pure Einstein-Hilbert truncation. The difference between the corresponding average 
values is approximately 1 for both 9' and 9" . 

Let us now come to the new critical exponent 9 3 which was not present in the Einstein- 



Hilbert truncation. Using the exponential shape functions ( 3.11 ) it suffers from relatively 
strong variations as the shape parameter s is changed. It assumes values in the range 
8.4 < #3(5) < 28.8. As compared to the exponential cutoffs, the cutoffs with compact 
support lead to a much weaker scheme dependence. For b 6 [0, 1.2] we have 23.0 < 9^{b) < 
26.7. However, the results obtained with the two families of shape functions agree within 
the scheme dependence. Moreover, #3 is always systematically larger than 9' (and 9") with 
both families of cutoffs. As a consequence, the hierarchy of critical exponents which was 
mentioned in (b) above and which squeezes the trajectories into a thin box is a universal 
feature. 

Obviously the critical exponents, in particular #3, exhibit a much stronger scheme de- 
pendence than g*A*. This is most probably due to neglecting further relevant operators 
in the truncation so that the B-matrix we are diagonalizing is still too small. 
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Figure 7: (a) 9' = Re#i and 9" = Im^i, and (b) 9% as functions of b, using the family of 
shape functions with compact support. 

In 2 + e dimensions 

The above results and their mutual consistency strongly suggest that 4-dimensional Quan- 
tum Einstein Gravity indeed possesses a RG fixed point with precisely the properties 
needed for its nonperturbative renormalizability or "asymptotic safety". However, with 
the present approach it is clearly not possible to determine the dimensionality Auv of 
the UV critical hypersurface, which coincides with the number of invariants relevant at 
the non-Gaussian fixed point. According to the canonical dimensional analysis, the (cur- 
vatures-invariants in 4 dimensions are classically marginal for n = 2 and irrelevant for 
n > 2. The results for 9% indicate that there are large nonclassical contributions so that 
there might be relevant operators perhaps even beyond n — 2. However, as it is hardly 
conceivable that the quantum effects change the signs of arbitrarily large (negative) clas- 



sical scaling dimensions, Auv should be finite [23 



A first confirmation of this picture comes from our .Recalculation in d = 2 + e where 
the dimensional count is shifted by two units. In this case we find indeed that the third 
scaling field is irrelevant for any cutoff employed, 9 3 < 0. 

For our analysis of the i? 2 -truncation in d = 2 + e dimensions with < e < 1 wc 
had to resort to numerical methods. Using the e-expansion we calculated the fixed point 
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coordinates and the critical exponents for selected values of the shape parameter s. For 
all quantities only the leading nontrivial order of the e-expansion was retained. In Table 
1 we present the corresponding numerical results. 



Table 1: Fixed point coordinates and critical exponents 


s 


X4+0(e 2 )) 


9* (+0(e 2 )) 


P* (+0(e)) 


9 X (+0(e)) 


02 (+0(e 2 )) 


0s (+0(e)) 


1 


-0.131e 


0.087^ 


-0.083 


2 


0.963^ 


-1.968 


5 


-0.055e 


0.092 £ 


-0.312 


2 


0.955^ 


-1.955 


10 


-0.035e 


0.095^ 


-0.592 


2 


0.955^ 


-1.956 



For all cutoffs used we obtain three real critical exponents, the first two are positive 
and the third is negative. Thus, the corresponding ^-direction is UV repulsive. This 
suggests that the dimensionality of S\jy could be as small as Auv = 2, but this is not 
a proof, of course. If so, the quantum theory would be characterized by only two free 
parameters, the renormalized Newton constant Go and the renormalized cosmological 
constant Ao, for instance. 

Let us now compare the results to those from the Einstein-Hilbert truncation [|]. |2|. 
The A- and ^-coordinates of the fixed point and the critical exponents 9\ and 82 are 
found to be similar to those in [I], ||. However, in the Einstein-Hilbert truncation the 
leading-order results g* = 3/38£: + C(£: 2 ) w 0.079 e + 0(e 2 ) and 9 2 = e + 0(e 2 ) are scheme 
independent, which is not quite true for the results above. Both truncations agree on 
9 1 = 2 + 0(e). 

Summary: 

Our main results concerning the non-Gaussian fixed point in the _R 2 -truncation are: 

(1^2) Position of the FP : The fixed point is found to exists for all cutoffs used. This 
result is highly nontrivial since the example of the Gaussian fixed point clearly shows that 
a fixed point of the Einstein-Hilbert truncation does not necessarily generalize to a fixed 
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point of the i? 2 -truncation. For every shape parameter the fixed point practically lies 
on the \-g— plane, and its position almost exactly coincides with that from the Einstein- 
Hilbert truncation. 

(2^2) Eigenvalues and -vectors: The fixed point is UV attractive in any of the three 
directions of the \-g-(3— space for all cutoffs employed. The linearized flow in its vicinity 
is always governed by a pair of complex conjugate critical exponents 61 = 6' + \6" = Q* 2 
with 9' > and a single real, positive critical exponent #3 > 0. It is essentially two- 
dimensional, and this two-dimensional flow is well descibed by the RG equations of the 
Einstein-Hilbert truncation. 

(3^) Scheme dependence: The scheme dependence of the critical exponents and of 
the product g* A* is of the same order of magnitude as in the case of the Einstein-Hilbert 
truncation. While the scheme dependence of 9" is weaker than in the case of the Einstein- 
Hilbert truncation we find that it is slightly larger for 9'. The exponent 0% shows a rela- 
tively strong dependence on the cutoff. The product g* A* again exhibits an impressively 
weak scheme dependence. 

(4^2) Dimensionality of iS>uv : The dimensionality Auv of the UV critical hypersurface 
cannot be determined within the present approach. However, the results from our R 2 - 
calculation in 2 + e dimensions suggest that Auv should be finite also in 4 dimensions. 

On the basis of the above results we believe that the non-Gaussian fixed point occuring 
in the Einstein-Hilbert truncation is very unlikely to be an artifact of this truncation but 
rather should be the projection of a fixed point in the exact theory. We demonstrated 
explicitly that the fixed point and all its qualitative properties are stable against the 
inclusion of a further invariant in the truncation. These results strongly support the 
hypothesis that 4-dimensional QEG is indeed nonperturbatively renormalizable. 



6 Positivity of action, Hessian, and cutoff 
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6.1 Positivity of the action 

It is a well known problem that in d > 2 dimensions the Euclidean Einstein-Hilbert action 

Seh[9) = ~ [ d d x^g{-R(g) + 2\} (6.1) 
167rG J 

is not bounded below. In fact, decomposing the metric as g^ u = exp(2x)g^ u where g^ u is 
a fixed reference metric we obtain 

S ^{9] = ~ J d d x^e^ [-R + 2Xe^-(d-l)(d-2)g^ (D» X ) (A*)] -(6.2) 

This shows that Seh can become arbitrarily negative if the conformal factor x{ x ) varies 
rapidly enough so that (D^x) 2 is large. Therefore it seems difficult to define a path 
integral Z = J T>g^ v exp(— S'eh) for Euclidean quantum gravity. 

The situation improves by including the term J d d x^/gR? with a positive coefficient 
since the resulting action is bounded below [18]. While the Einstein-Hilbert term J d d x^fg~R 
leads to a negative contribution to the kinetic term of the conformal factor, which domi- 
nates at small momenta, the i? 2 -term gives rise to a positive contribution dominating at 
large momenta. As a consequence, both the truncated action functional Tk[g,g] of eq. 
( [4.2| ) and the bare action S[g] = Ft^ [g,g] possess an absolute minimum. Moreover, 
rewriting the truncation ansatz (|4.2j ) with g^ u = g^ as 

-2\ 2 



T k [g, g] = Jd d x^g\ fa (r - ^) + W [aX, - 



(6.3) 



one can easily determine a sufficient condition for a manifestly positive action Tk[g, g] > 0. 
In terms of the dimensionless couplings it reads: > 0, [3k > 0, and 

1287r^A fc /3 fc >l. (6.4) 



6.2 Positivity of the Hessian 

At the level of the flow equation, Tk appears on the RHS in terms of its Hessian rj^ to 
which the cutoff operator TZk is adapted by the rule ( |3.5| ). Thus, only if 1^ is a positive 
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definite operator we can obtain a cutoff which leads to a "correct" mode suppression. 
Since we expect the i? 2 -truncation anyhow to be reliable only for large k, it is actually 
sufficient if rjj^ and 7Z k are positive definite for sufficiently large momenta p 2 = —D 2 . 
The reason is that, due to the factor dtIZ k (p 2 ) which emphasizes the region p 2 « k 2 , 
the traces on the RHS of the RG equation ( p.4| ) receive the dominant contributions from 
modes whose p 2 is close to k 2 . 

(2) 

In general F k [g,g] depends on both and the background metric g^ u . Here we 
concentrate on F k 2 \g,g] = F k with the two metrics identified. Furthermore, we assume 
that g^ u = g^ u is the metric of a rf-sphere with radius r since our projection technique 
requires these backgrounds only. In this case the eigenvalues p 2 = Ai(d, s) depend on the 
discrete quantum number I. The explicit expressions for Ai(d, s) are tabulated in appendix 
p. They are strictly monotonically increasing functions of I with lim^oo Ai(d, s) = oo. 

(2) 

In the following we show that the operator with k very large indeed becomes 
positive definite if it is restricted to the subspace spanned by the — D 2 -eigenfunctions 
with sufficiently large eigenvalues, certain assumptions on the couplings being made. The 
spherical harmonics T l ™, Tjf 1 , and T lm with I larger than a certain minimum value / m ; n 
provide a basis of this subspace. We shall concentrate on the conditions implied by the 
leading large-/ behavior. 

The Hessian F k [g, g] as given by the quadratic form ( |4.12| ) is a symmetric blockdiago- 



nal matrix. Therefore, according to the Jacobi criterion, the condition for positivity takes 
the s„„p.e form (if)^ > 0, (if )„ > 0, (if )„ > 0, and (if (lf>)„ - 

(rj^M > 0. For sufficiently large values of / the leading /-powers of Aj(d, s) are the 
dominating contributions to the entries of F k 2 ^ in eq. ( |4. 1 2| ) so that, in this limit, the 
above condition boils down to 

< (Ff)'^(Z Nk K 2 -p k R) A,(d,2) =}► Z Nk K 2 -p k R>0 (6.5) 

V / h T h T 



o < (ri 2) )__^2^AKf) =► ^!> (6.6) 
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< ( r i 2 %/^ 2 (^) 2 Pk(Ai(d,0)) 2 =► /3 fc >0 (6.7) 








^(A ; (d,0)) 3 >0 =► ZNkK " ^ > (6.8) 

a 

For nonnegative values of the gauge parameter, a > 0, this leads to the following restric- 
tions on the dimensionless couplings: 

g k >0, P k >0, k 2 /(32ng k [3 k ) > R. (6.9) 

In the UV fixed point regime of the d = 4-dimensional case we have g k ~ g* and 
/3 k ~ (3* with <7*,/3* > 0. Hence, close to the non-Gaussian fixed point, the first two 
conditions of eq. ( |6.9|) are obviously satisfied. Furthermore, the third condition then 
takes the form R < k 2 /(327rg*f3*). For R fixed this condition is satisfied as well provided 
k is sufficiently large. Thus, for k large and on modes with large eigenvalues of — D 2 , 

(2) 

the restricted operator T k is positive. The cutoff should have the desired suppression 
properties therefore. 

The above argument treats R as a constant parameter. Recalling the derivation of the 
projected flow equation where we compared powers of the radius r cx R~ x l 2 it is indeed 
clear that in this context r and R should be regarded as fixed, fc-independent quantities. 

It is instructive to look also at the operator rj^ [g os (k), g os (k)] where g os (k) is the 
fc-dependent "on-shell" S^-metric which solves the equation of motion STk/Sg^ = for 
9fiu = g^v The difference to the situation discussed before is that R is a function of k 
now, to be computed from g os (k). (The operator T k 2 \g os (k), g os {k)} would appear in a 
standard one-loop (saddle point) calculation based upon the "classical" action T k .) 

For the truncated action functional T k of eq. ( |4.2|) with g^ v = g^ u the field equation 
takes the form 

2Z Nk K 2 \Gftv + g^u Aj + h [~ (GV + R pv ) R + 2D^D U R - 2 9fMU D 2 R] = (6.10) 



d- 1 
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with G^u = R^v — the Einstein tensor. Precisely for d — 4, the maximally 

symmetric solutions to eq. ( |6.1U| ) satisfy Einstein's equation G^ u = —g^u^k, they are 
not affected by the i? 2 -term. Inserting the contracted equation R = AX k into the third 
condition of ( |6.9| ) leads to 

128ng k X k (3 k < 1 . (6.11) 

Remarkably, this condition is satisfied precisely if T k [g,g] is not a manifestly positive 
functional of g^ v , as follows from ( |6 . 4| ) . This implies that for d = 4 the S 4 solution of 
eq. fl6.10|) cannot correspond to the absolute minimum of T k [g,g] if this functional is 
manifestly positive. 

In the UV fixed point regime, the condition ( |6.11| ) becomes 1287T g*A*/3* < 1. For all 



cutoffs employed we found that 0.17 < 1287T g*\*/3* < 0.22 so that this condition is indeed 
satisfied. It is reassuring that also, upon inserting the S 14 solution of eq. fl6.10|) , the Hessian 

(2) 

T k becomes a positive operator for sufficiently large values of I and k, independently of 
the cutoff. 

Furthermore, the concomitant violation of ( |6.4| ) implies that in the vicinity of the fixed 
point the functional T k [g, g] is bounded below but not positive. By adding an appropriate 
constant it is trivial though to turn it into a manifestly positive functional. 



6.3 Positivity of the cutoff 

The cutoff A k S is expected to be positive definite under the same conditions as found 

(2) 

for the Hessian T k in the previous subsection. In order to obtain more quantitative 
information about Z m ; n and the momentum regime where A k S is positive we continue 
our analysis with an explicit investigation of the cutoff operator TZ k . For simplicity we 
again restrict our considerations to the most interesting case of d — 4 and to spherical 
backgrounds. 

After setting {U k )f T f T = 1/2 + g^g va ) (n k )- hT - hT and (K k )^ = <T (K k )& 

and inserting the eigenvalues of the covariant Laplacians, the entries of the cutoff matrix 



65 



1|) assume the form 

fc*fl(°)(A,(4,2)/fc 2 ) f in 
(**W = 32^; (1-32.^-1 

fc 4 i^°)(A,(4,l)/A; 2 ) 



(ft*) 



" 167T5ffca 
(W*)„ = ^ (0) ^ 0)A2) {^*9 A (2A ; (4,0)/A; 2 + ^°)(AK4,0)/fc 2 )) + jj} 

(tt*)* = (^ i = ^ fc A; 4 |[A,(4,0)A 2 + J R( )(AK4,0)A 2 )] 1 



x 1 /A,(4,0)A 2 + i2(0)(A,(4,0)/fc2) - A 



- [A ; (4,0)A 2 ] f ^(4,0) A 2 - Aj 1 j (6.12) 

faO* = ^ (0) ^1° )A2) {36vr^.(2AK4,0)A 2 + ^)(AK4,0)A 2 )) -1} . 

As compared to eq. ( |3.8|) which was written in terms of the dimensionful quantities Z^kK 2 , 
Afcj and we switched here to a description in terms of the dimensionless couplings. 

(2) 

In analogy with the Hessian T k the condition for the cutoff matrix 1Z k to be positive 

definite reads: (K k ) hT - hT > 0, (K k )g > 0, (K k )^ > 0, and {K k ) u (K k ) 9S - (TL k )\- > 0. 

(2) 

These conditions indeed reproduce the restrictions on the couplings obtained from T k 
for sufficiently large momenta p 2 = Ai(d,s). Provided that a > 0, they take the form 
g k > 0, (3 k > and k 2 / (327T g k f3 k ) > R, which coincides with eq. ( |6.9|) . 

Given an arbitrary set of parameters (R,k,g k ,/3 k ) satisfying these three inequalities, 
we have (^jO^t^t > and (TZk)^ > (and also (JZ k )-- > 0) for any allowed value of 
I. This is not the case for the other two conditions which stem from the scalar sector of 
the cutoff. Clearly (7£&)^ of eq. ( p.l2|) can assume negative values for sufficiently small 
values of I, provided g k , (3 k and R/k 2 are small enough. Since (JZ k )-- > 0, a negative 
(K k )iz implies that also (KJtt (JZ k )-- -(TZ k )%<0. 
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The /-values for which (Hk)$$ > satisfy 

2A,(4,0)/^ + fi<»)(A,(4,0)/F)> I5i i^. (6.13) 

In appendix |0J we derive a similar inequality involving the <p-a cross term. There we find 
that (7lk)$$ {Ttk)aa ~ {^k)\a > at least for all values of / satisfying 

2A,(4, 0) A 2 + fl<°> (a,(4, 0)/k 2 ) > — J— . (6.14) 
V / 96n g k (j k 

Both inequalities, ( |6.13| ) and ( |6.14j ), depend on g k [3 k , on p 2 /£; 2 with p 2 = A;(4, 0) = 
/(/ + 3)i?/12, and on the shape function R>°\ Given a specific set (R, k, g k , Pk] R^', T lm ) 
with a certain scalar eigenmode T lm , they tell us whether the contribution from this mode 
is suppressed correctly or not. The more restrictive inequality ( |6 . 1 4| ) applies to the scalar 
eigenmodes {T lm } with I > 2, while ()Sl3D can be used for the constant mode T ' m 
and the PCKV's {T l=1 > m } only. 

Let us now focus on RG trajectories which run into the non-Gaussian fixed point as 
k — > oo. Furthermore, we assume that R is either kept fixed or that R = a k k 2 with a 
constant a k < (327T gkPk)^ 1 - Then, for large enough values of k, we have g k ~ g* > 0, 
/3 k w /3„ > and A: 2 / (327T > so that the conditions ( |6.9| ) for the positivity of 

r[ 2) and are satisfied. Moreover, the RHS of (|6.13|) and ( |6.14D may be expressed as 
(1447T f?*/?*) -1 and (96n g^P^)" 1 , respectively. 

Now we are in a position to determine the /-regime for which the cutoff is manifestly 
positive definite. Using the family of exponential shape functions ( |3.il|) with 1 < s < 30, 
a numerical analysis reveals that any value of the ratio x = p 2 /k 2 = Aj(4, 0)/k 2 satisfies 
( |6.13j ) or ( |6.14| ) provided s > 2 or s > 7, respectively. Hence, under the above conditions, 
all cutoff's employing an exponential shape function with s > 7 are manifestly positive 
definite for all momenta, i.e. for all quantum numbers /. 

This is a rather intriguing result. It might indicate that cutoffs with s > 7 are 
particularly reliable. 

Conversely, for any s < 7 there exists a specific value Xo(s) such that all x with 
x < xq(s) violate ( |6.14| ). Furthermore, there exists a specific xi(s) for any s < 2 such 
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that any x < Xi(s) leads to a violation of ( |6.13| ). In FIG. § we show xq(s) and xi(s) in the 
ranges 1 < s < 6 and 1 < s < 2, respectively. It is important to note that Xq(s) < 0.7 and 
Xi(s) < 0.26 for any value of s considered. This implies that in the UV fixed point regime 
the cutoff has the desired suppression properties for all modes with momenta ranging 
from infinity down to values well below k. 
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(a) 

Figure 8: (a) x , and (b) x\ as functions of s, using the family of exponential shape 
functions. 



To complete the analysis let us study the inequalities ( |6.13| ) and ( 6.1 4j) also at the 
spherically symmetric stationary point g os (k) of Tk[g, g] which we discussed in the previous 
subsection. In the vicinity of the fixed point the on-shell value of the curvature is R ~ 
4A*fc 2 . Hence, we obtain from ( |6.13| ) and ( |6.14| ), respectively, 

1 



/l[jR (o). z] = 2l (l + 3)A, + R(0){1{1 + 3)K/3) 



144vr g^fc 



> , 1 = 0,1 



(6.15) 



and 



(0) ^ 2Z(Z + 3)A, + R{0) + 3 _ ]_ 

3 96ng*p* 



> , Z > 2 



(6.16) 



The first inequality stems from the scalar eigenmodes T lm with I = 0, 1, and the second 
from those with I > 2. Both Q6.15 ) and (|6.16|) depend on / and R^. 

Again we restrict our investigation to the family of exponential shape functions with 
1 < s < 30. Then the LHS of flOS]) and ( |6TT6|) are functions of s and /: /i[-R (0) ;Z] = 
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fi{s,l), f [R {0) ;l} = fo(sj). For I > 2 we have f {s;l) > f {s;l = 2) independently of 
the shape parameter. Numerically we find that fo{s; I = 2) is always positive. Hence, 
any momentum with / > 2 satisfies the condition fl6.16|) for all values of s considered. 
Furthermore, our numerical analysis shows that also fi(s,l = 1) > for all cutoffs 
employed. However, f±(s,l = 0) is not always positive. We obtain fi(s,l = 0) > for 
s > 2 and fi(s, I = 0) < for s < 2. Our results are illustrated in FIG. || 




(a) (b) 

Figure 9: (a) fo(s, I = 2), and (b) fi(s, I = 0) and fi(s, I = 1) as functions of s, using the 
family of exponential shape functions. 

The above results have to be interpreted as follows. Assume that k lies in the UV 
scaling regime, and consider the cutoff operator at the spherically symmetric stationary 
point, lZk[g os {k)]. Then this operator is strictly positive on the space spanned by all 
spherical harmonics and would correctly suppress all modes in a path integral containing 
this "on-shell" cutoff, provided we choose an exponential shape function with s > 2. For 
s < 2, only the contributions from the constant mode with p 2 = are not suppressed 
correctly. For any other mode the cutoff term is positive even in this case. 

To summarize: In this section we found that at least in the asymptotic domain 
relevant in our investigation of the UV fixed point the cutoff which is adapted to the R 2 - 
truncation is positive definite and therefore has all the required mode suppression proper- 
ties. No conformal factor problem and no growing exponentials produced by exp(— A^S) 
are encountered. 
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7 Summary and Conclusion 



In this paper we evaluated the exact RG equation of Quantum Einstein Gravity in a trun- 
cation which generalizes the Einstein-Hilbert approximation used so far by the inclusion 
of a higher-derivative term. We derived the beta-functions of the resulting \-g-/3— system 
which turned out to be by far more complicated than the old X-g— system. We used 
these beta-functions in order to investigate how the two fixed points known to exist in 
the Einstein-Hilbert truncation manifest themselves in the enlarged theory space. 

We found that the Gaussian fixed point of the Einstein-Hilbert truncation does not 
generalize to a corresponding fixed point of the i? 2 -truncation. Nevertheless, the 2- 
dimensional projection of the \-g-f3— flow onto the X-g— plane at (5 = 0, near the origin 
A = g = /3 — 0, is well approximated by the flow resulting from the Einstein-Hilbert 
truncation. The projected flow does indeed have a fixed point at A* = g* = 0. In the 
Einstein-Hilbert truncation there exists a distinguished RG trajectory, the "separatrix" 
f|, which gives rise to a vanishing renormalized cosmological constant, lirn^o Xk = 0. 
In d — 4, it turned out that this trajectory possesses a 3-dimensional "lift" which is 
characterized by a logarithmic running of the i? 2 -coupling For d ^ 4 its running is 
power-like, and there exists a "quasi- Gaussian" fixed point at A* = g* = 0, /3* ^ 0. This 
picture puts the older perturbative calculations in i? 2 -gravity into a broader context. 

Quite differently, the non-Gaussian fixed point (A*, g*) ^ (0, 0) implied by the Einstein- 
Hilbert truncation does "lift" to a corresponding fixed point of the 3-dimensional flow, 
with a tiny but nonzero third component /3* > 0. It is UV attractive in all 3 direc- 
tions of \-g-f3— space. We demonstrated in detail that close to the fixed point the flow 
on the extended theory space is essentially 2- dimensional, and that the 2-dimensional 
projected flow is very well approximated by the Einstein-Hilbert flow. For the A*- and 
(^-coordinates both truncations yield virtually identical values, and the same is true for 
the critical exponents pertaining to the 2-dimensional subspace. For universal quantities 
the differences between the two truncations are typically smaller than their weak residual 
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scheme dependence. 

This stability of the Einstein-Hilbert truncation against the inclusion of a further 
invariant, together with the other pieces of evidence which we summarized in Sections |5. 3. 1| 



and |5.3.2| strongly support the conjecture that this approximation is at least qualitatively 
reliable in the UV. Hence it appears increasingly unlikely that the very existence of the 
non-Gaussian fixed point is an artifact of the truncation. We believe that QEG has indeed 
very good chances of being nonperturbatively renormalizable. 

A notorious difficulty of Euclidean quantum gravity is the conformal factor problem. 
In the exact RG approach, it appears in the Einstein-Hilbert truncation, but not in the 
i? 2 -truncation provided k is large enough. When this complication occurs the construction 
of an appropriate cutoff operator is rather subtle. However, it was possible to show that 
our investigation of the non-Gaussian fixed point in the i? 2 -truncation is not affected 
by this problem, and that a straightforward positive definite cutoff can be employed. 
The numerical agreement of the results with those from the Einstein-Hilbert truncation 
indicates that the rule for constructing an adapted cutoff in presence of the conformal 
factor problem which was proposed in |l] ("Z^ = ^-rule") should indeed be correct. 
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A Evaluating the RHS of the truncated flow equation 



In this section we present several rather lengthy calculations needed for the discussion 
of the i? 2 -truncation in section |j. In the following, all calculations are performed with 
9fiu = 9fiu where g^ u is assumed to correspond to a S d background. For simplicity the 
bars are omitted from the metric, the curvature and the operators. 



A.l Computation of (vf ] + Kk) * and (s^ + Kk) 



In section f| we derived explicit expressions for the kinetic operators = V£ J + IZk and 



,(2) 



= + TZf.. They may be represented as matrix differential operators acting on the 
column vectors (h T , £, O , 0i) T and (v T , v T , g, g) T , respectively. In this representation 
they take the form 



?(2) 
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(A2) 



where 



Qk 



( 
\ 



V[g,g\ 



ri 2) [g,g] 
r k 2) [g,g} 



\ 
J 



(A3) 



The entries of these matrices are given in eq. Q4.15 ). On the RHS of the flow equation 

h 



the inverse operators [T^ ] 1 and [S^] 1 appear which we determine in the following. At 
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this point it is important to note that, because of the maximally symmetric background, 
all covariant derivatives contained in the operators ( |A1| ) and (|A^) appear as covariant 
Laplacians, and that the various entries are x-independent otherwise. This implies that 
these entries are commuting differential operators which allows for particularly simple 
manipulations. Therefore it is not difficult to verify that the inverse operators assume the 
form 



and 



r'i 2) [<?,<?] 



/ 
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Inserting these expressions into the RHS of the flow equation (|3.4j) leads to 



S k (R) = -Tr> 



1 . 







-1 


E 


(rf\g,g]+n k [g]) a 


dt(K k [g)) cc 







aa \ / 0101 V / 0icr 



;(2), 



^ \ / (TIT V / 0101 
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&1GT 
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where we used the relations 



s£W] + rc fc [(/] 
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(A7) 
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The trace of the 0o _ term appearing in the first line of eq. ( |A7| ) may be evaluated easily 

since only the scalar eigenmodes T 01 and T lm contribute. We obtain 

-i 



2 
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d t {n k [g\) 



1 Di(dfl) 



2Z ^ EE / ddx Vg&)T lm (x) [c S2 (d, a) C S i(d, a) (p k + A si (d, a) R 



1=0 m=l 
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+A S i(d, a)R + B S i(d, a) X k 

+ {Z Nk K 2 )- l f3 k (H s {d) (A l (d,0) + k 2 R(°\A l (d,0)/k 2 )y 
+G S2 (d) R (A,(d, 0) + k?RP\At(d, 0)/k 2 )) + G S3 (d) R 2 
xd t [c S2 (d, a) C sl (d, a) Z Nk K 2 k^R^^d, 0)/k 2 ) 
+P k (H s (d) {2A l (d,0)k 2 R^(A l (d,0)/k 2 ) + k i R^(A l (d,0)/k 2 ) 2 ) 



+G S 2(d)Rk 2 R {0 \A l (d,0)/k 2 



(A9) 



Here A;(d, 0) is the eigenvalue with respect to — D 2 corresponding to T lm . Inserting also 
the remaining entries given in eq. ( [4.15 ) into eq. ( |A7| ) finally leads to eq. ( [4.16 ). 



A. 2 Heat kernel expansion and evaluation of the traces 



In this part of the appendix we expand S k (R) of eq. ( OB ) with respect to r and evaluate 
the traces appearing in the resulting equation ( Alp) below by applying the heat kernel 
expansion. Thereby we extract the contributions proportional to J d d x yjg, J d d x ^/gR 
and J d x ^R 2 ■ This puts us in a position to read off the RG equations for the three 
couplings. For technical convenience we restrict our considerations to the gauge a = 1. 
We start our evaluation of S k (R), eq. ( |4.16| ), by expanding it with respect to R oc 



! . Since we are only interested in the contributions proportional to J d d x ^Jg oc r d , 



J d d x JgR 



oc r 



d-2 



and J d d x y/gR 2 oc r d 4 , only terms of order r d , r d 2 and r d 4 are 



needed. This leads to 

S k (R) = Tr (2ST2) [A^Af] + Tr' (1T) [A^Af] - h x {d) Tr' ( ' 0) [A 2 l N] + a k Tr'( 0) [A?%] 
+a*Tr' ( ' 0) [{A l A 2 )- 1 PlN] - 2Tr (1T) [P k ~ W ] - 2Tr' (0) [P^Af,] 
+ { - a fc Tr (25T2) [A^Ti] + a fc Tr (2ST2) [AfP k Af] - A T {d) Tr (25T2) [A^Af] 
-A v (d, l)Tr{ 1T) [AfAf] + h 2 (d) a k Tr'( 0) [(AiA 2 )' 1 T 2 ] 
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Here we set 
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= P k — 2X k , 
A 2 = a k P 2 - -h(d) (P k - 2~X k ) 
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The quantities A/", A/q, T\ and 72 are defined as in eq. ( |4.17| ). Furthermore, 0( 



,<d-4\ 



means that terms oc r n with powers n < d — 4 are neglected. 

The terms in eq. ( |A10|) proportional to 5^,2 and 5^,4 arise from the last term in eq. 
( |4.16| ). Contrary to the other terms of eq. ( |4.16| ), its expansion does not contain d- 
dependent powers of r, but is of the form b 2m r~ 2m with {62m} a set of r-independent 

m=0 

coefficients. As for comparing powers of r, this has the following consequence. Since, 
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for all m > and d > 0, —2m = d — 4 or —2m = d — 2 are satisfied only if (m, d) G 
{(0,2), (1,2), (0,4)}, and since —2m = d cannot be satisfied at all, this term contributes 
to the evolution equation only in the two- and the four- dimensional case. Using eq. 



the pieces contributing, i.e. &2m=o r° in d = 2 and d — 4, and &2m=2 r~ 2 in d = 2, may 
be expressed in terms of the operators J d 2 x^/gR, J d 2 x^fgR 2 or J d A x^fgR 2 . This yields 
the terms in eq. ( [A 10] ) which are proportional to the 5's. 

As the next step we evaluate the r-expansion of the traces appearing in eq. ( |A10|) by 
applying the heat kernel expansion. In its original form it has often been used to compute 
traces of operators acting on unconstrained fields. For our purposes we need the heat 
kernel expansions for operators acting on constrained fields, i.e. fields satisfying 
appropriate transversality conditions. In ref. these expansions are derived in detail 
for Laplacians D 2 on S d backgrounds acting on symmetric transverse traceless tensors, 
on transverse vectors and on scalars, with the following results: 
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x J i _ I (is _ e )fl _ 5 f ^ + 6 ( g + i£ ) 2 i? 2 + C)(i2 3 ) }> . (A15) 



6 V ' 360d(d-i; 
Here the terms proportional to the 5's arise from the exclusion of unphysical modes of 
the type discussed in section |2|. 

Let us now consider an arbitrary function W(z) with a Fourier transform W(s). For 
such functions W, we may express the trace of the operator W{— D 2 ) that results from 
replacing the argument of W with — D 2 in terms of W(s): 

oo 

Ty[W(-D 2 )] =lim J dsW(s) Tr e- {is - e)D2 (A16) 

— oo 

We obtain the asymptotic expansion of Tr [W(— D 2 )} by inserting the heat kernel expansion 
for Tr[e-( is ~ E ) D2 ] into eq. ( \KWj ). For Lapl acians acting on the constrained fields considered 



here they read as follows: 
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Q d/2 -i[W] y d d x^R 
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Tr (0) [W(-D 2 )] = (An)- d / 2 {Q d/2 [W} J d d x (A19) 
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i r ^d 2 — 7(7 -i- fi r ~\ 

+- Q Q d /2-i[W] J d d x JgR + ^ _ + Q d/2 - 2 [W] J d d x yfgR 2 + C(r <c!_4 )| . 
Here the set of functionals QnfH 7 ] is defined as 

oo 

Q n [W] = lim / rfs(-is + e)" n H/(s) . (A20) 

e\0 y 

— oo 

By virtue of the Mellin transformation we may now reexpress Q n in terms of W so that 

oo 

Q n [W] = Jr 1 ^., I dzz n + l ~ x , i > -n , iGNU {0} arbitrary. (A21) 

o 

In particular we obtain Qo[W^] = W(0). Furthermore, if n > we may choose % — for 
simplicity. As can be seen by an appropriate integration by parts, Q n [W] does not depend 
on i. 

At this point it is necessary to discuss the case where isolated eigenvalues have to be 
excluded from Tr[W(— D 2 )}. As we showed in ||, such traces can be expressed as 

Tt'-'[W(-D 2 )] = Tr[W(—D 2 )] - D ^ S ) W ( A ^' s )) • ( A22 ) 

ie{h,...,in} 

Here n primes at Tr' '"' symbolize the exclusion of all eigenmodes T lm with I G {1%, . . . , l n }, 
and Ai(d, s) and Di(d, s) denote the corresponding eigenvalues of — D 2 and their degrees 
of degeneracy, respectively. Since Ai(d,s) oc R we may view W(Ai(d,s)) as a function 
of R. As outlined above, such a function contributes to the evolution of Afc, Z^k and 
only for d = 2 and for d = 4, with the contributions given by W(0) + W'(0) Aj(2, s) and 
W(0), respectively. Using the explicit expressions for Ai(d, s) a.iadDi(d, s) (see Table 2 in 
appendix |C|) and applying eq. ( |4.7| ) we therefore obtain for the traces relevant to the flow 
equation: 

Tt' (1t) [W(-D 2 )] = Tr (lT) [W(-D 2 )] - ^ |eW(0) J d 2 x ^R 

+3W'(0) Jd 2 x^R 2 )-^yW(0) f d 4 x R 2 + O (r< rf - 4 ) , 

(A23) 
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Tr(o) [W(—D 2 )] = Tt (0) [W(-D')] - & |gW(0) | d 2 * ^J2 

+3W'(0) jd 2 x^R 2 ^- W(0) jd 4 x^R 2 + (r< rf - 4 ) , (A24) 

Tr{ 0) [W(-D 2 )]=Tr [0) [W(-D 2 )} 

J S W{0) I d2x ^ R ~2^f W{0) J R 2 + O (r <d - 4 ) . (A25) 

Here W denotes the derivative with respect to the argument: W'(z) = dW(z)/dz with 
z = A t (d,s). 

The next step is to insert the expansions of the traces into Sk{R), eq. ( |A10|) , and to 
compare the coefficients of the operators J d d x^fg 1 J d d x^/~g~R and j d d x^/g~R 2 with those 
on the LHS, eq. ( |4.5|) . This leads to the following differential equations: 

d t (Z Nk \ k ) = (4 K 2 )- 1 (47r)- d / 2 {/i 8 (d) Q d/2 [A^Af] - h^d) Q d/2 [A^Af] 

+a k Q d/2 [A 2 l T 2 ] + a k Q d/2 [(AxAz)' 1 P%M] - 2dQ d/2 [P fc "W ] } , (A26) 

d t Z Nk = -(2/ t 2 )- 1 (47r)- d / 2 {/i 9 ( C /)g d/2 _ 1 [A^M] - hMQd/2-1 [A^M] 

+-a k Q d /2-i [A^Tz] + ^a fc g d/2 _i [(AiA*)- 1 P£M] + h w (d)Q d/2 ^ [P fc _ W ] 
-h u {d) a k Q d/2 [A^] + h n (d) a k Q d/2 [AfP k N] + h 12 {d) Q d/2 [AfN] 
+h 2 (d)a k Q d/2 [(AiA 2 )~ 1 %] +h(d)a k Q d/2 [(AxA^ 1 P k Af] 
+h(d)Q d/2 [{A 1 A 2 )- l M}+h^d)a k Q d/2 [A^ 1 ^] -h 2 (d)a 2 k Q d/2 [A\ x A 2 2 P 2 k P 2 \ 
-/i 3 (d) a 2 Q d/2 [A 2 2 P k T 2 ] - /i 4 (d) a fc Q d/2 [A 2 2 T 2 ] - h 2 (d) a\ Q d/2 [{A X A 2 )- 2 P£M] 
-h 3 (d) 4 Q d / 2 [A[ l A 2 2 PtM] - 2h 4 (d) a k Q d/2 [A^A 2 2 P^] (A27) 
+^(d) /13(d) a fc Q d/2 [A 2 2 P k Af] + h(d) h 4 (d) Q d/2 [A 2 2 M] + h n (d) Q d/2 [P k 2 M ] } , 

d t pk = (4/ t 2 )- 1 (47r)- d / 2 {/i 14 (d) Q d/2 ^ 2 [A^Af] - ht(d) h 15 (d) Q d/2 . 2 [A^Af] 
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+h 15 (d)a k Q d/2 _ 2 [A 2 X T 2 \ + h 15 (d) a k Q d/2 _ 2 [{A^)' 1 PlN] 

-h w {d) Q d/2 _ 2 [P fc -Wo] - h 17 (d) a k Q dj2 _ x [Al x %\ + h 17 (d) a k Q d/2 ^ [AfP k M] 

-h w (d)Q d/2 -i [A{ 2 N] + ^h 2 (d)a k Q d/2 ^ [(A^)' 1 ^ 

+\h 3 {d)a k Q d/2 ^ [{AiA^PkN] + ^/i 4 (d)Qd/2-i [GAi^rW] 

+^ 3 (d)afcQd/2-i [^2^1] -\h 2 {d)a 2 k Q d/2 ^ [A{ 1 A 2 2 P 2 k ? 2 \ 

~h 3 (d)alQ d/2 -i [A 2 2 P k T 2 ] - ^h A (d) a k Q d/2 -i [A 2 2 T 2 ] 

—hMalQ^ [(A.A.y'P^-hMalQ^-i [A^A 2 2 PlM] 

~h A {d)a k Q d/2 -x [A^A 2 2 PlN] + ^h 1 (d)h 19 (d)a k Q d/2 _ 1 [A 2 2 P k U] 

+h l (d)h 4 (d)Q d/2 _ 1 [A 2 2 Af]-h 20 (d)Q d/2 ^ [P k 2 Mo\-h ll {d)a 2 k Q d/2 [AT 2 P k %] 

+^h u (d) h 21 (d) a k Q d/2 [A^] + h 22 (d) a k Q d/2 [A^ 2 Af] 
+h 11 (d)a 2 k Q d/2 [AfP 2 k M] -h u (d) h 21 (d)a k Q d/2 [AfP k M] 
+h 23 (d)Q d/2 [AfN]-h,{d)a 2 k Q d/2 [(A^)- 1 ^] 
+h 6 {d)a k Q d/2 [(AxA^N] -j 2 a\Q dl2 [{AxA^ 1 P k T x ] 
+h 2 (d) h 3 (d)a k Q d/2 [(M) _1 Ti] -h 2 (d) 2 a 2 k Q d/2 [(A^)' 2 P k T 2 ] 
+h 2 (d)alQ d/2 [A^A 2 2 P 2 T 2 ] -2h 2 (d)h 3 (d)a 2 k Q d/2 [A^A 2 2 P k T 2 ] 

--/12(d) h 4 (d) a k Q d/2 [A^A 2 2 T 2 ] - h 6 (d) a 2 k Q d/2 [A 2 2 T 2 ] 

+ \h 2 {d)alQ d/2 [{A l A 2 y 2 P k A N]-2h 2 {d) h 3 (d)a 2 Q d/2 [{A^)' 2 P^N] 

-\h 2 {d) h A (d)a k Q d/2 [{A 1 A2)- 2 PlM]-h T {d)alQ i/2 [A^AfP^M] 

-3h 3 (d)h 4 (d)a k Q d/2 [A^A 2 2 P k M] - hi,{d) 2 Q d/2 [A^A 2 2 N] 
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+h l {d) h 6 (d) a k Q d/2 [A 2 2 M] - h 2 {d) h 3 {d) a\ Q d/2 [A^AfP 2 ^] 
-h 3 (d) 2 a\ Q d/2 [A 2 2 P k %\ - h 3 (d) /i 4 (d) a k Q d/2 [A 2 2 %\ 
+h 2 (d) 2 a 3 k Q d/2 [A^ 2 A 2 3 P k i T 2 ]+2h 2 (d) h 3 (d)a 3 Q d/2 [A^A 2 3 P 3 T 2 ] 
+2h 2 {d)h A {d)a 2 k Q d/2 [A^A 2 3 P^T 2 ] + h 3 (d) 2 alQ d/2 [A 2 3 P^T 2 ] 
+2h 3 (d) h A (d) a\ Q d/2 [A 2 3 P k T 2 ] + h 4 (d) 2 a k Q d/2 [A 2 3 T 2 ] 
+h 2 (d) 2 a 3 Q d/2 [{A l A 2 )- 3 PtN]+2h 2 {d) h 3 (d)a 3 Q d/2 [A?A 2 *P%N] 
+3h 2 (d)h 4 (d)a 2 k Q d/2 lAfA 2 3 P£M}+h 3 (d) 2 a 3 k Q d/2 [A^A^P^N] 
+Ah 3 (d)h A (d)alQ d/2 [A^ l A 2 3 P 3 N]+3h,{d) 2 a k Q d/2 [A^A 2 3 P 2 k M] 
-h^d) h 3 (d) 2 a 2 Q d/2 [A 2 3 P 2 M] - 2h 1 (d) /13(d) h A (d) a k Q d/2 [A 2 3 P k M] 
-/11(d) /14(d) 2 Qd/a [A 2 W]+h 24 (d)Q d/2 [P k 3 K] 

+S d ,2 



3 *^ + (*RW {0) + l ak e 



2 ak p k (k*-2\ k ) 



— rW(q) dtZm 

4 1 } Z Nk (k 2 -2X k ) 2P 



T \2 



Z Nk (k 2 - 2X k ) 



^ (0) '(o) 



+5 



'<L4 



"I d t [Z Nk k 2 



1 



^ [Zivfc k' 



18 Z Nk (k 2 - 2X k ) 4 Z Nk (3a k k 4 - (k 2 - 2A fc )) 



-rO fe A; 4 



0> t [ZATfc fc 



cuzvfcA; 5 



8 ak Z Nk (k 2 - 2X k ) (3a k fc 4 - (A; 2 - 2X k )) 4 (9a fc A; 4 - 2 (k 2 - 2X k )) 



8 ak p k (3a k k* - (k 2 - 2A fc )) + 8 G " ft (9a fe A; 4 - 2 (k 2 - 2X k ))\ 



9 



}■ 



(A28) 



Here the coefficients hi are functions of the dimensionality d. They are tabulated in (|D2|) 



of the appendix |D.2 . 

Now we introduce the cutoff-dependent generalized threshold functions 



^"(v^d) 



(-iy 

T(n + i) 



dyy 



,n+i 



.1 # 
97 



(y + tf (0) (y))' 



y) -2/ 



(A29) 
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x (y + R(°\y) + W )~ P ( 32™ (y + R^(y)) 2 - JL-^ (y + R^(y) + w) 



and 



(y + R m (y)) m (2j/ + iJ<%))V(!/)(A30) 







x (y + Rf»(y) + w)- p ( 32n v (y + R^(y)) 2 - (y + Rp\y) + w) 



In eqs. ( |A29| ) and ( |A30| ), i is a nonnegative integer which satisfies i > —n, but which is 
arbitrary otherwise. The functions \1/ and \1/ are independent of i which can be seen by 
an integration by parts. Again, we may set i = if n > 0. Furthermore, noting that 



^ n (w) = ^(v,w;d), & n (w) = *™,o(v,w,d) Vv W (A31) 

we recover the threshold functions and $^(u>) originally defined in |l| in the context 

of the pure Einstein-Hilbert truncation. 
Using the relations 

Q^A^A^P^T,} = fc 2 ( m+ "-^ +1 )v]/^ m (a fc A; 2 /(32vr),A fc /A; 2 ;rf) 

-i|fe(A0 k^ m+n -?-^™ m . fi (a k fc 2 /(32vr), \ k /k 2 ; d) 

Q n [A^ p A 2 q P r k n T 2 ] = 2fc 2 ( m+ ^-^vl/^ m+1 (a fc A ; 2 /(327r),A fc /fc 2 ;rf) 

-l Vf) (k) k^ m+n '^ 2 ^ m . A (a k k 2 /(32n),X k /k 2 ; d) 

Q n [A7A^ q P^M] = k 2 ^ +n - p ^ +1 ^ m (a k k 2 /(32n),X k /k 2 -,d) 

~rj N (k) k 2 ^ +n ^ +1 ^ (a k k 2 /(32n),X k /k 2 ; d) 
Q n [p- p Af ] = k 2 ^ p+1 ^ p n (0) (A32) 



the differential equations (|A26 ), ( A27|) and ( |A28| ) may be rewritten in terms of the thresh- 



old functions ^ q m (v, w; d) and ^'.^(u, w; d) instead of the Q n . 
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In order to make the integrals in eq. ( |A20|) convergent we have to demand that (y) 
decreases rapidly as y — > ±00. However, since from now on its form for y < does not 
play a role any more we identify R^°'(y) with its part for nonnegative arguments and 
assume that R^(y) is a smooth function defined only for y > and endowed with the 
properties stated in subsection EO. 



Next we introduce the dimensionless couplings X k , g k and (3 k of eqs. (}4.20|) - (|4.22|) . 
Inserting eq. ( |4.20| ) into d t (ZN k X k ) leads to the relation 



d t Xk = ~ (2 - VN(k)) A fc + 32irg k K 2 k d d t (Z Nk X k ) . 



(A33) 



Then, by using eq. ( |A26| ) , we obtain the differential equation (|4.23|) for the dimensionless 
cosmological constant. The corresponding differential equations for g k and j3 k may be 
determined as follows. Taking the scale derivative of eq. ( |4.21D and ( ^.22| ) leads to eqs. 
( [4.24] ) and (|4.25|) , respectively. For the anomalous dimensions rjN and r]p we obtain from 
eqs. ( |M7D and (|A"28"D , respectively, 

i] N = g k Bx(X k , g k , f3 k ; d) + rj N g k B 2 {X k , g k , P k ; d) +r]j3g k B 3 (X k ,g k ,p k ;d) (A34) 

and 

VP = -Pk 1C ^k,9k,Pk,d) -VNp k 1 C 2 (X k ,g k ,p k ;d) - rj^ p^ 1 C 3 (X k , g k , p k ; d) . (A35) 



The set of equations (|A34|) and (|A35|) may now be solved for the anomalous dimensions 
rjN and rj/3 in terms of g k , (3 k and d which eventually leads to the expressions (|4.26|) 
and (lOTQ . 



B Coefficient functions appearing in the /3-functions 

In the following we list the coefficient functions Ai, Bi, Ci, i = 1, 2, 3, which appear in the 
/3-functions ( |4.23| ), ( }4.24| ) and ( |4 . 2 5| ) . The coefficients hi(d) contained in these expressions 
are defined in subsection 1072], and the generalized threshold functions \1/ and \1/ were 
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introduced in eq. (|A29|) and eq. ( |A30| ), respectively. The other threshold functions, $ 
and are those already introduced in the context of the pure Einstein-Hilbert truncation 

Ai(\ k , g k , p k ; d) = -2A fe + 2(4tt) 1 -5 g h (h s (d) ^ /2 (-2A fc ) - 2d $* /2 (0) 
+647T g k p k ^j 2 . x {gkPk, -2A fe ; d) - h x (d)^ 2 , (g k p k , -2X k ; d) 
+32tt g k f3 k VfypigkPk, -2A fe ; d) J (Bl) 



A 2 (A fc ,c/ fc ,/3 fc ;d) = A fe - (47T) 1 3 ih 8 (d) $J /2 (-2A fc ) 

+/i 1 (rf)$2) 1 2;0;0 ((7 fc /? fc , -2A fc ; d) + 32irg k /3 k ^} 2 . 2 , {g k Pk, -2A fc ; d) ) f B2 ) 



M*k,9k,Pk-,<I) = -8(4tt) 2 -2 gl/3 k ^} 2A1 ig k p k ,-2X k ;d) (B3) 

£i(A fc , d) = 4(4tt)H {/19(d) $i /2 _ 1 (-2A fe ) + /i 10 (d) ^/^(O) 

32 1 

+ ^ 9kPk^°Jl 2 -i;i(9kPk, -2A fc ;d) - ~/ii(d) ^^-^(sA -2A fe ; d) 



6 
16 

+ Y n 9kPk^ i /2 -. x . 2 {9kPk, -2A fc ;d) - 32/i 1 i(d)7r0 fc /3 fc $J /2 (-2Aj fc ) 



+32h 11 (d)7t g k p k ^ 2 . x (g k P k , -2A fc ; d) + h X2 (d) ® 2 d/2 (-2\ k ) + /i 13 (d) $* /2 (0) 
+32h 19 (d)n g k (3 k ^J 2 . x {g k P k -, -2X k ; d) + h A (d) ^a/ 2 . (g k P k , -2X k ; d) 
+32/i 3 (d)vr g k /3 k ^} 2 . (g k p k , ~2X k ; d) - h 19 {d) (32vr g k p k f ^J 2 . 3 (g k Pk, ~2A fc ; d) 
-QAh A (d)ix g k p k ^jJ 2 . 2 {g k P k , -2A fc ; d) - 2h 3 (d)(327T g k p k ) 2 ^a/ 2 . 2 (9kPk, ~2A fc ; d) 
+32h 1 (d)h 19 (d)ic g k p k ^ jJ 2 . x {gkPk, -2\ k ; d) + h x {d)h±{d) ^J 2 . (9kPk, -2\ k ; d) 
-h 2 {d){32irg k p k ) 2 ^ 2 . A {g k Pk, -2A fc ; d)| (B4) 
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B 2 (\ k ,g k ,p k ;d) = -2(4 7 r) 1 -i |/i 9 (d) $^ /2 _ 1 (-2A fe ) - -h^d) $};. 1;0 . ( Jt A, -2A fe ; d) 
16 ~ ~ 

+ ^9kPk *i/2-l;2;o(fl'^' ~ 2A ^ rf ) + 32/l n (rf)vr ^fc *5/2;l;0 ' ~ 2Afc ' rf ) 

+/i 12 (d) $3 /2 (-2A fc ) + 32h 3 {d)it g k (3 k *$ 2 . lfi (gkPk, ~2A fc ; d) 

+/i 4 (d) ^d} 2 -o;o(9kPk, ~2A fc ; d) - h 3 (d)(32rr g k (3 k ) 2 ^ 2 . 3 . (gkPk, -2\ k ; d) 

-64:h 4 (d)ir g k (3 k ^d? 2 . 2 . (g k Pk, -2A fc ; d) + 32h 1 (d)h 3 (d)7r g k (3 k *2/2;i ; o(^*/ 3 *' -2A fc ; rf ) 

+h 1 {d)h 4 (d) ^J 2 . . Q (g k Pk, -2A fc ; d) - /l 2 (d)(327T^) 2 ^/2;4;o(^^' ~ 2A ^ rf ) } ( B5 ) 

d f 16 ~ 

B 3 (\ k ,g k ,(3 k ;d) = -2(47r) 1 "5 | y7r<? fe /3 fe f ^^.^(ftA, -2A fe ; d) (B6) 
-32h u {d)it g k f3 k $^ /2 (-2A fe ) + 32h 2 (d)n g k (3 k #i) 2 . 0;1 ~2A fc ; d) 
+32h 3 (d)ir g k (3 k ^°j} 2 . . (gkPk: -2A fe ; d) - h 2 (d)(32ir g k [3 k ) 2 ^^(gkPk, -2A fc ; d) 
-/ i3 (d)(327T( ?fe A) 2 ^/2;i;i(rt> -2A fc ; d) - 32h i {d)ng k (5 k ^ 2;0;1 (^, ~2\ k ; d)| 

Ci(A fc ,^,/3 fc ;d) = (47r)-3 |/i 14 (d) $^ /2 _ 2 (-2A fc ) -/i 16 (d) $i /2 _ 2 (0) 

+6Ah 15 (d)n g k (3 k ^d} 2 _ 2]1 (gA, -2\ k ; d) - /ii(d)/ii 5 (d)^^ 1 2 _ 2;0 (5f fe /3 fe , -2A fc ; d) 

+32h 15 (d)n g k p k ^ 2 _ 2 . 2 (g k p k , -2A fc ;d) - 32h 17 (d)ir g k (3 k $^ /2 _ 1 (-2A fc ) 

16 16 

+^h 3 (d)n g k f3 k *°) 2 _ 1;0 (5'fe/3fe, -2A fe ; d) + y h l9 (d)ir g k f3 k *d/ 2 _ 1;1 (fffc/3fc, ~2A fc ; d) 

+ i/i 4 (d) *i) 2 _ 1 . (5'fe/3fe, -2A fe ; d) + 32h 17 (d)n g k f3 k ^a/ 2 -i;i(9kPk, -2\ k ; d) 
-h 18 (d) $ 2 /2 _ 1 (-2A fc ) - Md) ^/ 2 -i(0) 
-^i9(^)(327T5( fe /? fc ) 2 ^i) 2 2 _i ; 3(fi'fc/3fc, -2A fc ; d) 

1"2 2 0"2 

— jh 4 {d)7rg k /3 k ^^ 2 _ 1 . 2 {g k /3 k ,-2X k ;d) - -h 3 (d) (32k g k (3 k ) ^a/z-^KQkPk, -2A fc ; d) 



87 



X 6 1 

+ — h 1 (d)h 19 (d)7T g k (3 k *2/2-l;l W*' ~ 2X ^ d ) + g^l(<0/l4(d) *2/2-l;o(^/ 3 fc> ~ 2A ^ rf ) 

-i/i 2 (rf)(327T(7 fc /3 fc ) 2 ^ 2 2 _ 1;4 (<7 fe Ab -2A fc ; d) - /i n (d)(327r ^ fe ) 2 ^)° 2;1 (^/3 fe , -2A fc ; d) 

+Z2h 75 {d)Ttg k p k & d/2 {-2X k ) + /i n (d)(327r ^/? fe ) 2 *J) 2 . 2 (^/3 fc , -2A fe ; d) 

-32hn(d)h 21 (d)n g k (3 k ^ 20 {gkPh, ~2X k ; d) + h 23 (d) $ 3 d/2 (-2\ k ) + h 24 (d) $ 3 /2 (0) 

-h 2 (d)(32ir g k p k ) 2 ^} 2]1 (gA, -2X k ; d) + 32h 2 (d)h 26 (d)rt g k (3 k * d %. (g k Pk, ~2X k ; d) 

+h 2 (d){327T g k faf * d %. 3 (g k P k , -2X k ; d) - h 27 (d)(32n g k fa) 2 * d %. 2 (g k P k , -2X k ; d) 

3 

+32h 4 {d)h 28 {d)ir g k (3 k ^ d 2 2 . tl (gkPk, -2X k ; d) - -h 4 (d) 2 ^ d 2 2 . (gkPk, -2X k ; d) 
-h 29 (d)(32ir g k (3 k ) 2 ^^(gA, -2X k ; d) - 32h 4 (d)h 26 (d)ir g k (3 k V d J 2 . (g k P k , -2X k ; d) 
+ h 2 {d)(327i g k p k f * 2 d %. 4 (g k Pk, -2X k ; d) 
+hi 2 {d)h 28 {d){32ix g k f3 k ) 2 * 2 d %. 3 (g k P k , -2X k ; d) 

-48h 2 (d)h 4 (d)ir g k (3 k * 2 d %. 2 (g k P k , -2X k ; d) + 2/ i3 (d) 2 (32vr^) 3 *° d f 2 . 3 (g k P k , -2X k ; d) 

-3/i 1 (d)/ i3 (d)/i3o(rf)(327T g k f3 k ) 2 *° d %. 2 (g k Pk, ~2X k ; d) 
64 

+—h 1 (d)h 4 (d)h 28 (d)n g k (3 k V°t* 2il (g k Pk, ~2X k ; d) - h^h^df ® d %. (g k P k , -2X k ; d) 
+3h 3 (d)h 30 (d)(32n g k f3 k f ^ d %. 4 (g k P k , -2X k ; d) 

~h 4 {d)h 28 {d){32ix g k f3 k ) 2 V d %. 3 (g k p k , ~2X k - d) + §Qh 4 {d) 2 ix g k (3 k V d %. 2 (g k p k , -2X k ; d) 
- -h 2 (d) h 28 (d) (32tt g k k ) 3 ^ 3 2;5 (g k p k , -2 A fe ; d) 
+3h 2 (d)h 4 (d)(327ig k p k ) 2 * 2 d %. 4 (g k P k , -2X k - d) 
+h 2 (d) 2 (32irg k {3 k ) 3 * d %. 6 (g k Pk, ~2X k ; d) 

+ [2 (0) "T^2A7 + 2(l-2A fc ) 2 \ d ' 2 
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+ 



288tt^A-1 I 



4(1447r^/3 fc -(l-2A fc )) 4(l-2A fc ) 

96tt g k /3 k -l 24% g k f3 k 



2(96tt g k (3 k - (1 - 2A fc )) (1 - 2X k ) (96tt g k (3 k - (1 - 2A fc ) 



(B7) 



C 2 (X k ,g k ,p k ;d) = -^(4tt) 2 j/i 14 (d) $i /2 _ 2 (-2A fc ) - /i 1 (d)/ii 5 (d)*J 1 2 _ 2 . . (^/3 fc , -2A fc ;d) 

~ 16 ~ 

+32/i 15 (rf)7r5( A ./3 A; ^ d %_ 2 . 2 . (9kPk, ~ 2A ^ rf ) + y hz{d)-K g k (3 k ^^-i-i^k^k, ~2X k ; d) 

+^h A (d) *i/ 2 _i. 0i o(fffc/ 3 fc> _2 A fc ; rf ) + 32h 17 {d)-K g k {3 k ^a^i-i-oidA, ~2X k ; d) 

-h 18 (d) $^ /2 _ 1 (-2A fe ) - i/i 3 (rf)(327r^) 2 ^ 2 2 _ 1 . 3;0 (^/3 fc , -2A fc ; d) 

32 ~ 
— ^h 4 (d)ng k (3 k ^2-i;2;o(9kPk, -2A fe ; d) 



16 

+-^-/ii(d)/ii9(d)T^fc/5fc* d / 2 _i;i;o(^fc/ 3 fc> -2A fe ;d) 

1 



,<, 1V u,,,<, 19W „ y kHk -* d/2 _ 



+-h 1 (d)h 4 (d) ^° d %_ 1 . , (gkf3 k , ~2X k ; d) - -h 2 (d)(32n g k {3 k ) 2 ^/2-i;4;o(^fe, ~2A fc ; d) 



+32/i22 (d)ng k (3 k ® 2 d/2 (-2X k ) + h u (d)(327r g k p k ) 2 * 3 d %. 2 . (g k (3 k , -2A fc ; d) 
-32h n (d)h 21 (d)7vg k p k ^.^(gA, -2A fc ; d) + h 23 (d) $d/ 2 (-2A fc ) 

+32h e (d)ir g k p k ^ d } 2 . . (gkPk, -2X k ; d) - h 7 (d)(32ir g k (3 k ) 2 ^ d f 2 . 2 . (gkPk, ~2X k ; d) 

~ 3 ~ 

-9Qh 3 (d)h A (d)ir g k [3 k *J) 2 2;1 . (fffc/3fc, ~2A fc ; d) - -h 4 (d) 2 ^ d %. . (g k P k , ~2A fc ; d) 

+32h 1 (d)h 6 (d)TT g k f3 k ^ d f 2 . . (gA, ~2X k ; d) 

1 

2 



+ i/i 2 (d)(327T^/3 fe ) 3 ttj 2 2;4 . o (0 fc /3 fc , -2A fc ; d) 



-2h 2 (d)h 3 (d)(32ir g k /3 k ) 2 ^ d %. 3 . (g k (3 k , -2X k ; d) 
-48/i 2 (d)/i 4 (d)7r g k (3 k ^ d J 2 . 2 . {gkPk, ~2X k ; d) 



-h 1 (d)h 3 (d) 2 (32ng k (3 k ) 2 ^° d %. 2 . (g k (3 k , ~2X k ; d) 



(2/2;: 
0:3 



-64/ii(d)/i 3 (d)/i 4 (d)7r g k (3 k ^ J /2 . ilfi (9kPk, -2A fe ; d) 
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-h^h^d) 2 *% 2m (g k Pk, -2A fc ; d) 
+h 3 (d) 2 (32irg k [3 k ) 3 ^/ 2 . A . (g k p k , ~2X k ; d) 
+4h 3 (d)h 4 (d)(32irg k (3 k ) 2 ^ 2 . 3 . (gA, -2\ k ; d) 
+9Qh 4 (d) 2 ir g k (3 k ^ 2 . 2 . (g k p k , -2X k ; d) 
+2h 2 (d)h 3 (d)(32ir g k p k f *% 2 . 5 . (g k Pk, ~2\ k ; d) 
+3h 2 (d)h 4 (d)(32n g k (3 k ) 2 ^L. {gA, -2X k ; d) 



+h 2 (d) 2 (32irg k p k ) 3 * 3 j/ 2 . 6 . (g k P k , -2X k ; d) 



+ 



+ 



UR(°y (0) l_92ng k (3 k + 11^°)' (0) 



2(l-2A fc ) 



1 



2(l-2A fc ) ; 



+ 



J (2,2 



1 



4(U47rg k (3 k -(l-2\ k )) 4(1 - 2\ k ) 

1 24irg k (3 k 



2(96ng k (3 k - (1 - 2A fc )) (1 - 2A fc ) (?&K 9k p k - (1 - 2A fc )) 



(B8) 



C 3 (\ k ,g k ,p k ;d) = -i(4vr) a |32/z 15 (c/)7r 5f fc /? fc ^ 2 _ 2 . . 1 {g k P k , -2A fc ; d) 

-32h 17 (d)ir g k (3 k $^ /2 _ 1 (-2A fc ) + —h 3 (d)n g k f3 k *^2-i;o ; o(ff*/ 3 *' _2Afc ' rf ) 

16 ~ 
+ ^-h 2 (d)irg k [3 k ^a/2-i-,o;i(9kPk, -2A fc ; d) 



-^ 2 (rf)(327r^ fe ) 2 V^^gkPk, "2A fc ; d) 



-^ 3 (rf)(327rrt) 2 ^S) 2 2 _i ; i;i(^A, -2A fc ; d) 



16 



— ^-h 4 (d)7r g k p k yj /2 _ 1 . . 1 (g k Pk, -2A fc ; d) - /in(d)(327r^ fc /3 fc ) * d ) 2 ;i;o(^*' _2Afc ' rf ) 
+16^n(d)/i2i(d)7r^ fc ^ /2 (-2A fc ) - h b {d){32rx g k p k f ^l^kPk, -2A fc ; d) 



^(327r^/5 fc ) 2 *% 2il . (gkPk, -2A fc ; d) + 32h 2 (d)h 3 (d)7r ^/? fc * 1 j / 1 2 . . (gkPk, ~2A fc ; d) 



i;i 
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+ -h 2 {d){32ngA) 3 ^ 2 . 2;1 (gA,-2X k ;d) 



-h 6 (d) (32n g k p k ) 2 ^a^.^igA, _2 ^! rf ) _ h 3 (d) 2 (32n gA) 2 ^° d f 2 .x ;0 {gA, ~2X k ; d) 



-2h 2 (d)h 3 (d)(32ngA) 2 ^.^(gA, ~2X k ; d) 
-48h 2 (d)h 4 (d)iTgA ^^.^(gA, ~2X k ; d) 
-h 2 (d)h 3 (d)(32n gA) 2 ^ d %. 2 . (gA, ~2X k ; d) 

0;2 
d/2;l 

-32h 3 (d)h A (d)ng k p k ^ 2 2 . Q . (g k (3 k , -2X k ;d) 

-h 2 (d) 2 (32ng k {3 k ) 2 & d %. 2 . x {gA, -2X k ; d) + h 3 (df ; (32tt gA) 3 ^.^(gA, -2X k ; d) 
+2h 3 {d)h A {d){32ix g k f3 k ) 2 ^.^{gA, ~2X k ; d) 
+32h 4 (d) 2 ir gA ^^{gA, ~2X k ; d) 
+2h 2 (d)h 3 (d)(32ngAf ^.^(gA, -2X kl d) 
+2h 2 (d)h A (d)(327r gA) 2 ^ 2 . 2 . tl (gA, -2X k - d) 

2;3 
d/2;. 



+/i 2 ((i) 2 (327r^ fc ) 3 *d/2- 41 (^A, ~2X k - d) 



96TtgA 
' 1 - 2X k 



367r ^ (3 k 



2Att g k f3 k 



lUng k (3 k 



2X k 



96ng k p k 



2X k 



Jd,4 



(B9) 



In eqs. (|B7|), (|B5| ) and (|B^) the terms proportional to 5 dj2 or 5 d 4 arise not only from 



the 5-terms of eq. (|A10|) , but also by evaluating the "primed" traces, i.e. by subtracting 
the contributions coming from unphysical modes, see appendix |A.2| for details. All these 
contributions are obtained by expanding various functions f(R) with respect to R and 
retaining only the terms f(0) + f'(0)R in d = 2 and f(0) in d — 4. As we explained above, 
these are the only pieces of / which may contribute to the evolution in the truncated 
parameter space. Furthermore, the heat kernel expansions of the traces corresponding to 
differentially constrained fields introduce additional contributions proportional to 5 dt2 or 
5 dA into eqs. (©-(ll)- 
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C Tensor spherical harmonics on S 

The spherical harmonics Tj^, Tjf 1 and T lm for symmetric transverse traceless (ST 2 ) ten- 
sors h^ v , transverse (T) vectors and scalars <p on S d form complete sets of orthogonal 
eigenfunctions with respect to the covariant Laplacians. They satisfy 

-D 2 T^(x) = At{d,2)Tp(x) , 
-D 2 Tl m {x) = A^,l)^ m (x), 

-D 2 T lm {x) = A l {d,0)T lm {x) (CI) 
and, after proper normalization, 

5 lk 5 mn = j 'd d xV§ (l(2ST^r Pa T l ™T% = j d d x^g (l (1T) )^ T l ™ T kn 

= J d d x^/g~T lm T kn . (C2) 

Here (lpsT*)) 1 """ = (d - 2)/(2d) + 9^9 vp ) and (lfiT))" = (d-l)/dg^ are 

the unit matrices in the spaces of ST 2 tensors and transverse vectors, respectively. The 
Ai(d, s)'s denote the eigenvalues of —D 2 where s is the spin of the field under consideration 
and I takes the values s, s + 1, s + 2, ■ ■ -. The index m = 1, • • • , Di(d, s) is a degeneracy 
index. 

In ref. El explicit expressions for Aj(d, s) and the degeneracies Di(d, s) were derived 
which are summarized in Table 2. The eigenvalues are expressed in terms of the curvature 
scalar R = d(d — l)/r 2 of the sphere with radius r. 

The spherical harmonics TJ£, T l ™ and T lm span the spaces of ST 2 tensors, T vectors, 
and scalars so that we may expand arbitrary functions h/F v , £ M and according to 

oo D t (d,2) 

h U*) = E E h T™T l ™(x), 

1=2 m=l 
oo D t (d,l) 

u*) = E E fc™^*). 

Z=l m=l 
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oo Di(d,0) 

= EE <PimT lm (x). (C3) 

1=0 m=l 

Eqs. (|C3|) may now be used to expand also any symmetric non-T 2 tensor and non- 
transverse vector in terms of spherical harmonics since they may be expressed in terms 
of ST 2 tensors, T vectors and scalars by using the decompositions (|2.3j ), ( p.7|) , see e.g. 

I!§ n, 0, g. 

Note that the £>i(d,l) = d(d + l)/2 modes {T^" 1 } and the £>i(d,0) = d+1 modes 
{T 1,m } satisfy the Killing equation ( |2.5|) and the scalar equation (|2.6|) , respectively, and 
that T ' 1 = const. Arbitrary symmetric rank-2 tensors receive no contribution from these 
modes. In the case of arbitrary vectors the constant scalar mode does not contribute. 
Such modes have no physical meaning and have to be omitted therefore. 



Table 2: Eigenvalues of — D 2 and their degeneracies on the ci-sphere 


Eigenfunction 


Spin s 


Eigenvalue Ai(d, s) 


Degeneracy Di(d, s) 


1 


T l ™{x) 


2 


l(l+d-l)-2 p 


(d+l)(d-2)(l+d)(l-l)(2l+d-l)(l+d-3)\ 


2,3,-.. 


d(d-i) n 


2(d-l)!(/+l)! 


T l ™(x) 


1 


i(i+d-i)-i p 
d(d-i) rL 


l(l+d-l)(2l+d-l)(l+d-3)\ 

{d-2)\(i+iy. 


1,2,-.. 


T lm (x) 





l(l+d-l) p 
d(d-l) rL 


(2l+d-l)(l+d-2)\ 
l\(d-l)\ 


o,i,... 



D Tables of coefficient functions 

(2) 

D.l Coefficients introduced in 1^ [g,g] 

In this subsection we define the various A's, B's, C's and G's and Hs{d) which appear in 
eqs. ( |4.12fl - (|4.16|) of subsection fO| and in eqs. (|A9|) and ( |A10| ) of appendix [A| 



. (J . d{d-3)+A d(d-5) + 8 . . a(d-2)-l 

d(rf-l) ' 2rf(rf-l) ' Md ' a) = d 

i j\ d - 4 , . aid - 4) 

G v (d) = — , A sl (d,a) = 



2d ' ° v ' ' 2a(d - 1) - (d - 2) 
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a i i \ a (d - 2) - 2 , 2ad 

^ S2 (d, a) = — ^— — J — — - , B si (d, a) - 



a(d-2) -2(d-l) ' ° v ' y 2a(d - 1) - (d - 2) ' 

r (A \- 2ad r (A \- 2a(d-l)-(d-2) d-2 

552Ka) = «(d-2)-2(d-l) ' ^ lKa) = ~4(d-l)-2a(d-2)d-l 

,_d-1 2(d-l)-a(d-2) (d-2)(a-l) 
652(0, a) = —5 , C 53 (d, a) 



d? a ' °° v ' ' a(d-2) -2(d-l) ' 

(d-l)(d-4) (d-l)(d-6) 
^si(d) = ^2 ' G 52W = j 2 , 

(d_-4)(d-6) „^_ 9 /rf-lV mn 
D.2 Coefficients appearing in the /3-functions 

Next we define the coefficients hi(d) contained in the /3-functions (|4.23|) , Q4.24|) and (|4.25|) 
via the coefficient functions Ai, Bi, Ci, i = 1, 2, 3, given in appendix [B|. They also appear 
in the approximate solutions for the non-Gaussian fixed point of appendix [Ej. 

h(A\- d - 2 h (A\- d ~ A h(A\- d ljzM±± h ( A\ — (d-2)(d-4) 

h{d) = d^T ' h ^ = — • = 2d{d _ x) > Hd) = — , 

, , „ d 2 - 4d - 2 / 7\ (d - 4) 2 / 7\ 5d 4 - 48d 3 + 148d 2 - 112d + 16 
h{d) = 2d 2 ' W " 2d(d^T) ' h7id) = 4d 2 (d-l) 2 ' 

_ d 2 + d-4 _ (d + 3)(d + 2)(d 2 -5d + 2) d 2 - 6 

w = — 2 — ' h9{d) = i2d(d-i) ' hw{d) = ~^r> 

, . (d + l)(d-2) , . d 4 - 2d 3 - 5d 2 + 16d - 14 

Mrf) " — 2 — ' Mc/) " htt) ' 

, , n 2(d + 1) , , . 5d 6 - 7d 5 - 139d 4 - 545d 3 - 898d 2 + 504d - 360 
hi 3 {d) = , hu{d) - 



2 



d ' " v 7 720d 2 (d- 1) 

_ 5d 2 - Id + 6 _ 5d 4 - 7d 3 - 54d 2 - 180d + 180 

Md)= 360d(d-l) ' hw{d)= 180d 2 (d-l) ' 

(d + 2)(d+l)(d-5) (d + 2)(d 5 -5d 4 -5d 3 + 43d 2 -68d+ 18) 

hn(d) = — — — — , li 18 {d) 



12(d- 1) ' lov ' 12d 2 (d- l) 2 

. , _ 5d 2 - 28d + 20 (d + 3)(d-2) 
Md)= 2d(d-l) ' Md) = 3d 2 " ' 
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{d)=2 d(d-l) ' Md) = 4d(d^lj ' 

d 6 -5d 5 + 3d 4 + 31d 3 -86d 2 + 98d-50 d + 3 

i(rf) = 2d 2 (d-l) 2 ' Md) = "2-^-, 

> {d) = AdJd^T) ' Md)= d(d-i) ' 

15d 4 -178d 3 + 628d 2 - 632d + 176 9(d 2 -6d + 4) 

' (rf)= 4^(^-1)2 > Md)= 2d{d _ l) , 

, A 5d 4 - 52d 3 + 168d 2 - 128d +16 , . A 3d 2 - 16d + 12 
i(o) = jirn — Tv2 ' h m d ) = 



id) 

i(d) 
,(d) 
t (d) 
5 (d) 



t (d) 



4d 2 (d-l) 2 ' " uv 7 2d(d-l) 

5d 6 - 27d 5 - 71d 4 - 405d 3 - 342d 2 - 960d + 360 
720d 2 (d- l) 2 ' 

d 6 - 3d 5 - 7d 4 + 5d 3 + 26d 2 - 82d + 12 
12d 2 (d- l) 2 ' 

d 6 - 5d 5 + 7d 4 - 13d 3 + 42d 2 - 42d + 2 
2d 2 (d-l) 2 ' 

5d 6 - 7d 5 - 119d 4 - 593d 3 - 846d 2 + 480d - 360 
360d 2 (d- l) 2 ' 

d 6 - 3d 5 - lid 4 + 9d 3 + 54d 2 - 134d + 36 



,(d) = 3- 



3d 2 (d- l) 2 
5d 5 + 7d 4 - 9d 3 + 46d 2 - 62d + 14 



2 



d 2 (d- 1) 

/ n _ (d + 2)(d 3 -6d 2 + 3d-6) d 4 -2d 3 + 3d 2 -Ad -2 

' (rf)= 3d(d^l) ' h ^ d > = - 2 dJdZT) ' 

. 5d 2 - 7d + 6 , . . 30d 5 - 115d 4 - 362d 3 + 721d 2 + 182d + 264 

'(") = Arjfj — ?rr~ ' h M d ) 



45d(d-2) ' 4UV 7 90d(d- l)(d-2) 

,3d 6 - 17d 5 + 25d 4 + 39d 3 - 166d 2 + 224d - 96 
' 3d 2 (d-l)(d-2) ' 



(d-l)(d-4) 2 {d + 2){d 3 -6d 2 + 3d-6) 
l[d)=A d{d-2) ' HaM= 3d(d-l)(d-2) 

/A d 4 -2d 3 + 3d 2 -4d-2 , , „ d 4 - 3d 3 + 32d - 32 
l( ' } " 2 d(d-l)(d-2) - kM = d 2 
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, /A d 4 - 13d 2 - 2Ad + 12 , . d 4 - 2d 3 - d 2 - Ad + 2 

6dJd^T) ' Mrf) " Kd^T) • (D2) 



E Closed-form formulas for the fixed point location 

In the following we derive the approximate formula for the position of the non-Gaussian 
fixed point discussed in subsection |57|. Here we restrict our considerations to the case 
d > 2. 

In a first approximation we set X k = A* = 0, (3 k = (3* = and determine g* from the 
condition //at* = 2 — d alone. Since /3* = 0, we may solve this equation for g* in closed 
form which leads to 

9 * ~ B 1 (0, 0, A*; d)-{d- 2)5 2 (0, 0, A*; I) ' (E1) 
As A* = 0, it boils down to 



g t = (Arc)- 1 \h A3 (d) $Vi(°) + /l 46(rf)^/ 2 -i(0) + h M (d) & d/2 (0) 



+Md)*2/ 2 m • (E2) 



Here the hi(d) are again <i-dependent coefficients which are defined in subsection p.2| of 
appendix |B|. It is remarkable that the solution ( |E2| ) coincides precisely with the corre- 
sponding approximate solution (H2) of ref. || with a = 1, obtained in the framework of 
the Einstein-Hilbert truncation. 



Employing the exponential shape function ( |3.11 ) with s = 1, and setting d = A 



for instance, eq. ( [E2|) yields g* ~ 0.590. Here we used that, for this shape function, 
$i(0) = tt 2 /6, $2(o) = l, $J(0) = 1, $2(0) = 1/2, see appendix 0. 

In order to improve upon this approximation scheme, we determine (A*, from 
a set of Taylor-expanded /3-functions. Using eqs. (|FlJ) -(|F4|) we expand the /3-functions 
O), and (|425l) about A fe = g k = (3 k = and obtain 

/3 x (\ k ,g k ,i3 k ;d) = -2\ k + v d d g k + O (g 2 ) , 
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(3 g (\ k ,g k ;a,d) = (d - 2) g k - (d - 2)u d g 2 k + O (g 3 ) 
^(0,0, /3 k ;d) = ld + (4-d)f3 k + 0{g 2 ) . 



(E3) 



Here 7^, v dy and Ud are defined as in eqs. fl5.8|), ( |5.12| ) and ( |5.18| ), respectively, and (9 (g n ) 
stands for terms of nth and higher orders in the couplings gi(fc) = A^, g2{k) = g k and 
g 3 (/c) = (3 k . Now g* is obtained as the nontrivial solution to j3 g = 0, which reads 



9* 



wj 1 = (4VT)!- 1 {h i3 (d) ^ /2 _ x (0) + M<0 <^ /2 (0)} 



-1 



(E4) 



Inserting eq. (|E4|) into (3 X = leads to 
d d(d — 3) j 



A* 



$} /2 (o) {/143(d) ^ /2 -i(o) + ^(d) $ jU(o)} 



-1 



(E5) 



Quite remarkably, also these results agree completely with those of ref. which follow 
from the pure Einstein-Hilbert truncation. (See eqs. (H6) and (H7) of this reference.) 

Now we use f3* in order to determine /3*. However, since the term linear in (5 k vanishes 
for d = 4, the expanded (3p of eq. ( |E3| ) is not sufficient in this case. Therefore we consider 
also those terms of second order in the couplings which are linear in (3 k . For these terms 
we find 



d 2 (3 



Od = 



d\ k d(3 k 



d 2 (3 fl 



dg k d(3 k 



Afc=9fe=/3fc=0 



Afe=9fe=/3fe=0 



0. 



-(4 7 r) 1 -f {2h 39 (d) + 2h A5 (d)<S> 2 /2 (0)-3 5 dA } • (E6) 



Taking the nonvanishing term of eq. (|E6| ) into account, and inserting g* of eq. (|E4j) into 
(3p = then leads to 

Id 



0: 



d- 4 - adU, 1 



(E7) 



(47T)-! \ 


^i/a-aCO) + /132(d) $3/2-1(0) + /133(d) 3^(0)' 


i 


d-4+ < 


2h 39 (d)+2h 45 (d) $ 2 d/2 (0)-3 5 dA 


M 


"^43(d)^ /2 _ 1 (0) + /144(d) ^ /2 (0)} _1 
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Employing the shape function ( |3.11|) with s = 1 we obtain from eqs. (E4|) , (E5|) and 
( E7| ) in d = 4 dimensions 

H£ + IlV" M 0.751 , 
144 24tt / 

= 419(13,' + 474) 
M (4vr) 2 906768 v ; 

Here we used the expressions for the threshold functions derived in appendix |F|. The 
numbers in ( |E8| ) should be compared to the exact result ( |5.33| ). 



F Properties of the threshold functions 

In this appendix we summarize various important properties of the threshold functions 
*n™> *n-L> $n and $p which are defined by eqs. flggg ), (p0|) and flX3| ). 



Expanding the generalized threshold functions ^1^, ^'m-z about vanishing couplings 
yields 



x ^+«-m + i (0) Afc + 2 ( p + q )( P + q + 1) (-2^) " $^-™+ 2 (0) A 2 fc 
-327rg [-2— ) ^ n +<1 ~ m ~\0) g k (3 k + 0(g 3 ) , (Fl) 



x^ m+1 (0) X k + 2(p + q )(p + q + 1) (-2^)" $r^ m+2 (0) X 2 k 

(i -, \ q+l 
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+2( P +g) (-2^) q ($r 9 - m (o) +n$ p „:r m+1 (o)) a, 

+2( P + q )(p+ q +i) (-2^) 9 ($r^ m+i (o) + n $ p „:r m+2 (o)) a* 

-327rg ("2^) ^ (^ +9 — 2 (0) + n ^ST^CO)) ifcft + ^(g 3 ) • (F3) 

Here (9 (g 3 ) stands for terms of third and higher orders in the couplings gi(fc) = A&, g2(&) = 
gk and g3(fc) = /3fc. Quite remarkably, every fixed order of these expansions depends only 
on the "conventional" threshold functions and at vanishing arguments. 

By using eq. (|A31|) the corresponding expansions of $^ and $^ about vanishing 
argument can be read off directly from eqs. ( |FT|) and ([F2j). They are given by 

$£(-2A fc ) = <S>l(0) + 2p^\0)\ k + 2p(p+l)^ 2 (0)\l + O(\l) 
K(-2X k ) = $l(0) + 2p^\0)\ k + 2p(p+l)& n + 2 (0)\l + O(\l) (F4) 

For n = the threshold functions are universal in the sense that they do not depend 
on R {0) (y). In fact, setting n = in ty™ m , *^.„ and leads to 

w; d) = w; d) = (l + w)- p (s2ir v - (1 + w)) ' (F5) 

and 

= = (1 + w)~ p . (F6) 

There exists a second class of universal values of certain threshold functions. Using 
the boundary conditions for R^(y) one may easily verify that, for vanishing argument 
and for n + 1 = p > 1, assumes the universal value 

•WO) = • (FT) 
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Let us now be more specific and opt for the family of exponential cutoffs ( |3.11| ). In this 



case the integral that defines $^ can be carried out analytically for vanishing argument. 
Using the integral representation of the polylogarithm 



1 C x z 71 ^ 1 

U n (x) = — dz- z , (F8) 

1 (n) J e z — x 



o 



one obtains for p — 0, ... ,4: 

$°(0) =n(n + l)s- n ((n + l) (F9) 

$^(0) = ns~ n {((n + 1) - Li n+1 (l - s)} (F10) 



*n(0) ={ ^ ) U ) (Fll) 

1 s = 1 



$3(0) = (F12) 

' [2(n-\)(\-sfY x s*- n {(2-s)Li n _ 2 (l-s)-sLi n _ 3 (l- S )} n^l, s^l 
[2 n - 1 (n-l)]~ 1 (2 n - 1 -l) "^1, s=l 

-[2(1 - s) 2 }- 1 s 3 - n |(2-s)LiL 1 l 0) (l-s)-sLiL 1 f ) (l-s)} n = l, s^l 
ln(2) n = s = 1 

$4(0) = (F13) 

[6(n - l)(n - 2)(1 - s) 3 ] -1 s 4 " n {2(s 2 - 3s + 3) Li n _ 3 (l - s) 

-3s(2 - s) Li n _ 4 (l - s) + s 2 Li n _ 5 (l - s)} n ± 1, 2 , s ± 1 

[(n - l)(n - 2)] _1 (1 - 2 3 -" + 3 2 -™) n ^ 1, 2 , s = 1 
[6(2n - 3)(1 - s) 3 ]" 1 s 4 -™ |2(s 2 - 3s + 3) Li^Cl - s) 

-3s(2-s)Li^ 4 ) (l-s) + s 2 Lii 1 i ) (l-s)} nG {1,2}, s^l 

ln(27/16) n = s = 1 

k ln(4/3) n = 2 , s = 1 
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Here we defined 

Li^)(z) ee ^/yLi„(x) (F14) 
an K ax 1 

and used the relations 

Li n (l) = CW, U^(x) = l ±^ (F15) 

x 

with C denoting the Riemann zeta-function. For nonvanishing arguments an analytic 
solution to the integrals defining the threshold functions is not known. 

For the exponential cutoff ( |3.11| ) with s = 1 there even exists a very useful relation 



among $^(0) and $^(0). One may easily verify that 

^(0) = $r 1 (0). (F16) 
This relation allows us to calculate the <3>^(0)-integrals analytically as well. 

G Proof of the inequality ( |6.14| ) 

In this appendix we prove the inequality ( |6.14| ). As a first step we consider the function 



/(a, y) = a^a (a - y) - y/l-y (Gl) 

with a > 1 and < y < 2/5. (For a = 1 this function vanishes identically: f(l,y) = 0.) 
An upper bound for this function may be obtained as follows. For the first two derivatives 
of / with respect to y we obtain 

= #-/(a,y) = -— =^==+ 1 



dy Jy 2^a (a-y) 2^1^' 

f i0 ' 2) (a,y) ee ^-J(a,y) = w + L_ . ( G 2) 

d V 4 [a (a-y)] 5 4(l-y) 5 

Solving f( 0,1 \a,y) = for y leads to the single solution 

_ a(a + 1) 

y = yo= 2 z n ■ G3 

cr + a + 1 
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Since f (0 > 2 \a,y ) = (a 2 + a + l) 5 (a 3 - 1) /(4a 3 ) > 0, we have f(a,y ) < f(a,y) for all 
y G [0, 1] and a > 1. Hence, for a > 1 fixed but arbitrary, / monotonically decreases in 
the interval y G [0,yo] where yo < 1. 

Furthermore, yo — 2/o(o) is a monotonically increasing function of a for all a > 1. 
Therefore we have that ?/o( a ) > 2/o(o = 1) — 2/3. As a consequence, / monotonically 
decreases in the interval < y < 2/3 for any value of a > 1. Since we restricted our 
considerations to y G [0,2/5] we obtain f(a,y) < f(a,0). 
(TZ k )j- may be obtained from / by replacing 

A ; (4,0)A 2 + ^°)(A,(4,0)A 2 ) r 4 

AK4,0)A 2 ' ^3A,(4,0) Z(Z + 3)' 1 } 

and multiplying the result by 9/3 k k 4 (Aj(4, 0)/4/c 2 ) 3,/2 . Note that for all I > 2 we have 
A;(4, 0) > 5R/6 so that y < 2/5 is indeed satisfied. Moreover, a > 1 is satisfied as 
long as J R(°)(A i (4,0)//c 2 ) > 0. If i2 (0) (A*(4, 0)/k 2 ) = 0, the cutoff in the scalar sector 
is zero anyway: (Tl k )^ = (JZ k )^- = (Tl k )^ = 0. Hence, for positive values of (3 k , 
/(a, y) < /(a, 0) leads to 

(W*)*»< ^A; 4 J R( )(A / (4,0)A 2 ){2A / (4,0)A 2 + J R(°)(A i (4,0)A 2 )} . (G5) 

Next we insert (7£ fc )-- and {TZ k )^ of eq. ( |6.12| ) into (TZ k )^ (TZ k )-- - {TZ k )^- and then 
use ( |G5|) . This yields 



(^»)# - > ^ — 32^r /,w j (G6) 

x |l87r^ fc (2A / (4,0)A 2 + i?W(A i (4,0)A 2 )) - A 

Obviously the positivity condition (TZk)^ (R- k )-- — (JZk) 2 ^- > now boils down to 

v(k 2 ,l,R) + 2AK4,0)A 2 + J R(°)(A i (4,0)A 2 ) > Qfi 1 - . (G7) 

Here v is the nonnegative function of k, I and R which represents the contributions to 
- {KkYfa neglected on the RHS of Q. We see that eq. ( |6.14|) is a sufficient condition 
for the inequality ( |G7| ) to be valid, i.e. for (7£fc)x£ {Jl k )-- — (TZ k ) 2 ^- to be positive. This 
is what we wanted to proof. 



k 4 R^(A l (4,0)/k 
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